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Chapitre 1 

Introduction 



Les hypergraphes, une generalisation des graphes [9] , sont des structures discretes, 
lis permettent une description et un niveau d'abstraction necessaires pour la concep- 
tion et l'analyse en algorithmique [26^ [30] . Nous notons, dans ces references, le recours 
aux series generatrices pour faire de l'analyse enumerative. Nous souhaitons obtenir des 
caracteristiques quantitatives sur les hypergraphes, pour cela nous utilisons pour l'essen- 
tiel les series generatrices qui permettent l'enumeration exacte et asymptotique selon la 
taille de ces structures. Karohski M. et Luczak T. , dans [2D], etudient les hypergraphes, 
nous nous distinguons de leur travail par l'utilisation des series generatrices offrant une 
concision aux preuves. 

Le plan de la these est le suivant : 

- Dans le second chapitre, nous parlons d'enumeration exacte des hypergraphes conn- 
exes. 

- Dans un troisieme chapitre, nous procedons a l'enumeration asymptotique de ces 
structures. 

- Et dans le quatrieme chapitre, nous etablissons quelques caracteristiques des hy- 
pergraphes deduites de la performance de Palgorithme glouton d'hypercouplage ou 
deduites du processus d'hypergraphe evoluant. 
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Chapitre 2 

Enumeration exacte 



Dans ce chapitre, nous adoptons deux manieres d'enumerer des composantes classees 
selon Vexces une relation liant le nombre d'hyperaretes et le nombre de sommets. Une 
premiere maniere est remuneration bijective : en mettant en evidence une bijection entre 
les structures. Une seconde maniere de proceder a remuneration des hypergraphes est 
celle que nous qualifions de recursive. Une telle distinction est aussi adoptee par Wright 
E.M dans [22] et [23] pour les graphes, lui permettant d'un cote de justifier la forme des 
series generatrices et de l'autre d'automatiser le calcul de ces series. Avant de proceder a 
ces enumerations, nous precisons quelques definitions et notions. 

2.1 Definitions et notions 

Definition 2.1.1. Un hypergraphe est un couple (V,S) , un ensemble V de sommets et 
un ensemble £ d'hyperaretes soit de sous ensemble de V . 

La plupart du temps, sauf mention contraire, un hypergraphe dans cette these est 
b-uniforme, c'est a dire que chacune de ses hyperaretes contient b sommets. Ainsi, la 
structure d'hypergraphe generalise la structure de graphe qui est alors vue comme un 
hypergraphe 2-uniforme. 

Definition 2.1.2. L'exces d'un hypergraphe H = (V,£) est 

exces(W) = ^(| ei |-l)-|U|. (2.1) 

L'exces d'un hypergraphe TL = (V, £) (fc-unifbrme) est 

exces(W) = (b - 1)\£\ - \V\ . (2.2) 

Cette notion d'exces, utilisee dans [21] , permet de classer les hypergraphes. Par exemple, 
nous avons les definitions suivantes : 

Definition 2.1.3. Un hyperarbre est une composante ou hypergraphe connexe d'exces 
— 1 , valeur minimum. 
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Definition 2.1.4. Un hypercycle est une composante d'exces 0. 

Definition 2.1.5. Une composante est dite complexe si elle est d'exces £ > 1 . 

2.2 Enumeration bijective 

Dans ce chapitre, pour enumerer les hypergraphes, nous choisissons de nous focaliser 
aux structures connexes et de distinguer les structures selon leur exces, nous procedons 
alors a l'enumeration des structures connexes des plus "simples" aux plus "complexes" 
dans le sens ou les plus simples sont les hyperarbres d'exces —1 , viennent ensuite les 
hypercycles d'exces puis les composantes complexes d'exces I > 1 donne dans l'ordre 
croissant de ce dernier. Dans cette section, nous adoptons un point de vue bijectif en 
exhibant clairement une bijection ou en explicitant a travers les series generatrices une 
telle bijection. 

2.2.1 Vue bijective des hyperarbres 

Dans le cas des graphes, il y a plusieurs manieres d'enumerer les arbres (voir [U [29j 
[T2"l [2j [16]) en particulier via le code de Priifer que nous generalisons ici afin d'enumerer 
les forets d'hyperarbres (ou d'arbres) enracines. 

Definition 2.2.1. Une foret d'hyperarbres enracines est un ensemble non ordonne d'hy- 
perarbres enracines. 

Definition 2.2.2. Une feuille est un groupe de (b — 1) sommets (non racine dans le cas 
de structure marquee) de degre 1 dans une meme hyperarete. 

La connaissance du nombre des forets d'hyperarbres enracines est un outil cle pour 
enumerer les structures qui peuvent etre decrites de maniere "concise", c'est a dire que 
les structures sont simplifiees en elaguant recursivement les feuilles. En particulier, nous 
serons amenes a considerer les structures ainsi elaguees selon leur taille qui sera le nombre 
d'hyperarbres enracines contenus dans la foret definie par l'elagage. 

Pour enumerer les forets a (k + 1) hyperarbres enracines, ayant n sommets et s hy- 
peraretes, nous procedons a leur codage comme un quadruplet (R, r, P, N) ou 

R : un ensemble de (A; + 1) sommets parmi {1, . . . , n} 
r : un sommet de R 
< P : un partitionnement non ordonne de s sous-ensembles, (2-3) 

chacun de taille (6 — 1) , de {1, . . . , n}\R 
,7V: un (s - l)-uplet de {1, . . . , n} s_1 . 

Dans le code d'une foret d'hyperarbres enracines, nous pouvons facilement lire : 

- le nombre de composantes qui n'est autre que le nombre de racines soit \R\ , 

- un sommet racine r qui rattache la feuille de la derniere hyperarete dans le processus 
d'elagage , 

- le nombre |P| = dim(iV) + 1 d'hyperaretes , 
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Algorithme 1 : Codage d'une foret d'hyperarbres enracines. 
Entrees : Une foret, de (k + 1) hyperarbres enracines, ayant s hyperaretes et 

n = s(b — 1) + k + 1 sommets. 
Sorties : Codage (R, r, P, N) defini par ([23]) . 
debut 

(R,r,P,N) <- ({racine},r, {},()) 
repeter 

Ajouter l'ensemble des sommets de la plus petite (dans l'ordre 
alphabetique) feuille dans le partionnement P et placer le sommet qui le 
relie dans le tirage N . 

Redefinir la foret sans les sommets de la plus petite feuille. 
jusqu'a plus aucune hyperarete dans la foret 

(Le dernier sommet place dans le tirage est necessairement une racine). 
Definir r comme le dernier sommet placer dans le tirage. 
retourner (R, r, P, N) 

fin 
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(R,r,F,N) : 

' R = {5,9, 13, 16} 
r = 13 

( P = {{1, 22}, {2, 17}, {3, 19}, {4, 8}, {6, 7}, {10, 15}, {11, 18}, {12, 14}, {20, 21}} 
N = (21,18,13,13,4,18,21,7) . 

Fig. 2.1 — Une foret d'hyperarbres enracines et son code. 

- le nombre d'hyperarbres non reduit a leur racine correspondant au nombre des 
racines distinctes apparaissant dans (N, r) . 
Dans l'algorithme [TJ de codage, dans la boucle, parmi les feuilles, le choix de la plus pe- 
tite feuille dans l'ordre alphabetique nous permet de retrouver sans ambiguite la foret. 
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En effet, nous avons l'algorithme [2] de decodage qui retourne une foret car lorsque les 
hyperaretes sont iterativement formees, les hyperarbres existants se grandissent, se re- 
joignent et finissent par s'accrocher a une racine en gardant leur structure d'hyperarbre. 
Plus precisement, un sommet x de N accroche un ensemble de la partition de P avec la 
garantie que x soit plus proche d'une racine par rapport aux sommets, de Pensemble de 
la partition, qui necessairement finissent par etre connectes a une racine (eventuellement 
x est une racine). 



Algorithme 2 : Decodage en foret d'hyperarbres enracines. 
Entrees : Des entiers naturels n, k et s verifiant n = s(b — 1) + k + 1 codage 

(R,r,¥,N) defini par (1231) . 
Sorties : Foret d'hyperarbre enracine, dont les sommets racines sont les (k + 1) 
sommets de R . 

debut 

repeter 

Former une hyperarete avec le premier sommet du tirage ./V et avec les 
sommets du premier (l'ordre etant celui induit par les etiquettes) ensemble 
du partitionnement P ne contenant aucun sommet qui apparait encore dans 
le tirage restant. 

Supprimer du partitionnement l'ensemble utilise et supprimer de la liste du 
tirage le premier sommet. 
jusqu'a liste du tirage vide 

Former avec le dernier sous-ensemble du partitionnement et avec le sommet r 
le s-ieme hyperarete. 
retourner la foret obtenue 

fin 



Theoreme 2.2.3. Le nombre de forets de (k + 1) hyperarbres enracines ayant s hy- 
peraretes est : 

Hk + 1) ) /U ' n'- 1 , (2.4) 



jfc + y v J [(b-iy.] s s\ 

avec le nombre de sommets n = n(s) = s(b — 1) + k + 1 . 

Preuve. La preuve decoule directement de la bijection que definit l'algorithme [T]de codage 
de telles forets par les ensembles de quadruplets (R, r, P, N) defini en (|2,3p . 

En fixant k = dans ce theoreme, nous obtenons 

Corollaire 2.2.4. Le nombre d'hyperarbres enracines ayant s hyperaretes : 

[(6-l)!]- S ! n ' (2 ' 5) 
avec le nombre de sommets n = n(s) = s(b — 1) + 1 . 
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II est immediat que ce result at generalise le result at dans le cas des arbres en prenant 
6 = 2. 

Corollaire 2.2.5 (Cayley). Le nombre d'arbres enracines ayant n sommets est n n_1 . 

Une avantage que presente une telle demonstration bijective est la possibility d'effec- 
tuer une generation aleatoire afin d'apprendre quelques caracteristiques des structures 
etudiees. Nous laissons cela comme perspective, pour le moment. Grace au theoreme 
12.2.31 nous sommes aussi en mesure de donner une expression explicite du nombre des 
hyper cycles. 

2.2.2 Vue bijective des hypercycles 

Les hypercycles sont des composantes les plus simples apres les hyperarbres dans le 
sens oil ces structures sont d'exces . L'enumeration des graphes unicycles est associee 
a Alfred RENYI dans [3j . Dans le cas plus general des hypergraphes unicycles, pour 
proceder a l'enumeration de ces structures (la meme idee est applicable pour les com- 
posantes complexes a exces £ > 1 fixe), afin d'utiliser le theoreme 12.2.31 . les structures 
sont recursivement elaguees jusqu'a ce qu'il n'y ait plus aucune feuille et obtenir ainsi des 
structures lisses. Ainsi, a un hypercycle ayant n sommets, nous ferons correspondre de 
maniere naturelle un couple (B-,F) avec B une structure lisse de taille (k + 1) et T , une 
foret de (k + 1) hyperarbres enracines ayant n sommets, obtenue a partir du processus 
d'elagage. Les etiquettes de B seront canonisees dans {1, . . . , k + 1} en respectant l'ordre 
croissante des etiquettes. 




Fig. 2.2 - Hypercycles 3-uniformes lisses non etiquetes de longueur respectivement 2 , 3 
et 4. 

Remarque 2.2.6. Comme le processus d'elagage consiste a supprimer les sommets des 
feuilles ainsi que les hyperaretes qui les accrochaient, il produit une structure de meme 
exces. 

Nous retrouvons le resultat, associe a Selivanov dans [T7], pour l'enumeration des 
hypercycles : 
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Theoreme 2.2.7. Le nombre d'hypercycles ay ant s hyperaretes est 

n!n s_1 (6 — 1 



2[(6-l)!]« f^si(a-j) 
avec le nombre de sommets n = n(s) = s(b — 1) . 

Preuve. Un hypercycle de longueur de cycle j correspond a une foret de j(b — 1) hy- 
perarbres enracines a un arrangement pres de ces derniers pour les differentes fagons de 
former le cycle. La foret aurait (s — j) hyperaretes et j(b — 1) composantes a arranger en 
cycle. Le nombre d'hypercycles ayant une longueur j de cycle et ayant s hyperaretes est 
done 

n \ . 1} [(s-j)(b-l)]l fl [7-(6-l)]! ] (2?) 

j(b-l)) 3{b lj [(6-l)!]-i( s -i)! n /12[(6-2)!P7' [2J) 

avec n = n{s) = s(b— 1) : le premier facteur entre accolade denombre des forets d'hyper- 
arbres enracines et le second denombre les hypercycles lisses etiquetes avec {1, . . . , j(b—l)} 
ayant j hyperaretes. En simplifiant cette equation, nous trouvons le terme de la somma- 
tion du theoreme, et comme la longueur de cycle peut prendre toute valeur entre 2 et s , 
nous obtenons le resultat en sommant sur j . (} 

Ainsi, pour faire la preuve du theoreme nous etions amenes a distinguer les hypercycles 
selon la longueur du cycle. Une question a poser est : pour un nombre de sommets 
n = s(b — 1) , quelle est la longueur j de cycle de la classe (selon j) qui contribue le plus 
au nombre des hypercycles ? Nous laisserons cette question en perspective. 

Par rapport au resultat d'enumeration des hyperarbres (voir des forets) , l'expres- 
sion du nombre des hypercycles est beaucoup plus complexe car requiert la sommation. 
Neanmoins, l'expression de ce nombre reste explicite. Expliciter le nombre des compo- 
santes complexes parait etre une tache ardue. 

Dans la suite, nous precisons comment determiner le nombre de composantes, en 
decrivant la bijection via les series generatrices qui permettent une lecture directe des 
operations combinatoires sur des structures. 



2.2.3 Introduction aux series generatrices exponentielles 

A une sequence de nombre (a n ) n6 ]N , nous associons la serie generatrice exponentielle 
(SGE) 

oo n 
n=0 

Les SGEs servent pour remuneration de structures avec une etiquette propre a chaque 
sommet et l'indice n , dans Pecriture ci-dessus, definit la taille de la structure etiquetee. 
Une grande avantage de l'utilisation des SGEs est la facilite de lecture qu'elles offrent : 
a partir des operations sur les series, nous sommes en mesure d'interpreter en terme 
d'operations sur les structures et inversement. donne le dictionnaire suivant pour 
faire cette lecture : 
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\J\J£;l dljlUIlo oUI luo oil UtlUIco 


S( i.H, pnTVPcnnn ri q n't" 
UVJij I Co UUllLLcU.lL 


.Au S 


A(z) + BlV) 


AxB 


A(z) ■ B(z) 


Substituer dans A par B 


AoB(z) 


Sequence de A 


(l-A(z)y 1 

A(z) k 
k\ 


Groupe de k A 


Ensemble de A 


exp (A(z)) 


Cycle de A 


-ln(Vl-^)) 


Marquage de k sommets de A 


fc! H^ A i z ) 



Nous nous servirons aussi de Toper ateur sur les series : 



(2- 



avec A(z) defini a (12. 8p . Familiarisons nous avec l'utilisation des SGEs a travers les 
exemples de structures rencontrees jusqu'ici. 



Definition 2.2.8. La SGE T des hyperarbres enracines est 

m = E 



s>0 



[(6-l)!] s s! ra! 



(2.10) 



oil n = n(s) = s(6 — 1) + 1 . 

Nous y lisons le nombre des hyperarbres enracines a s hyperaretes 

(n- 1)! 



ra! [z n ] T(z) 



[(&- l)!] s s! 



(2.11) 



Remarque 2.2.9. La SGE qui enumere les hyperarbres s'ecrit sous la forme -£/_i oT(z) 
Nous reservons la justification pour plus tard. 

Un marquage d'un sommet differencie les hyperarbres des hyperarbres enracines : 



T(z) = z^-H^oT(z) 
Un cycle de sommets de longueur > 2 admet la SGE 



ln(Vl - z) - z/2. 



(2.12) 
(2.13) 



Pour un cycle de sommet avec un ensemble de (b — 2) , structure de SGE z b 1 /(b — 2)! , 
la substitution dans l'equation precedente donne 



yb-l 



J>-1 



In 



1 



(6 — 2)! 7 2(6-2)! 



(2.14) 



Et comme c'est la SGE des hypercycles lisses, pour obtenir la SGE des hypercycles, il 
faut faire la substitution par la SGE des hyperarbres enracines. De cette maniere, nous 
obtenons : 
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Remarque 2.2.10. Si b > 3 , la SGE des hypercycles est Hq o T(z) avec 



H (t) = -\n 




1 - 



(6-2)! 



) 



2(6-2)! ' 



(2.15) 



Nous pouvons lire en sens inverse la bijection qui nous avait permis de proceder a 
l'enumeration des hypercycles dans la fonction #0 • Cette fonction offre aussi la possibility 
de determiner le nombre n! [z n ] Hq o T(z) d'hypercycles de taille n . 

Comme il a ete deja evoque, il est possible d'appliquer, comme dans le cas des hyper- 
cycles, un raisonnement proche pour faire de l'enumeration des composantes complexes. 
C'est ce que nous allons developper dans la prochaine sous-section. 

2.2.4 Vue bijective des composantes complexes 

Dans cette sous-section, pour enumerer les composantes complexes, nous renongons a 
etablir une expression explicite comme dans le cas des forets d'hyperarbres enracines 
ou comme dans le cas des hypercycles. Nous y etablissons la SGE des composantes 
complexes d'exces donne, par des arguments bijectifs comme pour ([2.150 . Dans le cas 
de l'enumeration des composantes de graphe complexes, l'enumeration des composantes 
ayant deux cycles (d'exces 1 ) est associe a Bagaev dans [25] . Plus generalement, un 
resultat pour les composantes de graphes d'exces I > 1 est associe a Wright dans [22J . 

Ici, nous enumerons des composantes (d'hypergraphe) complexes d'exces £ > 1 . Pour 
proceder a cette enumeration, une premiere etape consiste a se ramener a des composantes 
lisses par elagage des feuilles des structures a enumerer. Les composantes lisses ainsi 
obtenues ont le meme exces que les structures de depart. 



Nous convenons d'utiliser la SGE avec la variable t pour enumerer des structures 
lisses, comme dans ([2.150 . Et dans le but d'alleger les equations, notons 




Fig. 2.3 - Une composante lisse non etiquetee d'exces I = 1 . 



V 



6-1 





(2.16) 
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(|2.15p devient alors 

H (t) = - In (y/T=7®) - ^ • (2.17) 

Les composantes complexes lisses d'exces I peuvent avoir toute taille positive (entiere) . 
C'est aussi le cas pour les hypercycles et nous avions eu recours a la fonction In du 
dictionnaire pour exprimer la structure cyclique. Dans le cas des composantes complexes, 
c'est l'operation de sequence du dictionnaire qui nous permettra de decrire les structures 
de maniere "plus concise" car pour un exces I donne nous aurons seulement un nombre 
fini de structures a considerer pour tout enumerer. 

Definition 2.2.11. Une chaine est une sequence de sommets avec un groupe de (b — 2) 
sommets. La sequence peut etre vide et admet pour SGE 

T^wy < 2 ' 18 > 

Une chaine, dans cette these, est d'abord une sous-structure d'une composante lisse. 
Dans le cas des graphes, Wright E.M dans [22] distingue trois types de chaine 

- celle qui boucle sur un meme sommet, 

- celle qui, si brisee, deconnecte la structure en deux, 

- celle qui ne boucle pas sur un meme sommet et, qui si brisee, ne rompt pas la 
connexite de la composante. 

Une chaine du premier (respectivement du second (respectivement du troisieme)) type 
pour des graphes sans multi-arete admet au moins deux aretes (respectivement une arete 
(respectivement deux aretes)). 



Definition 2.2.12. Une structure basique est une structure lisse non etiquetee ayant ses 
chaines de longueur finie. 

Rappelons ici l'idee cle dans [22] de l'enumeration bijective des graphes connexes 
d'exces i > 1 : a partir des composantes, une reduction est faite par elagage recursive, 
puis les chaines sont effondrees jusqu'a leur taille minimale et enfin, les etiquettes des 
sommets sont ignorees pour obtenir un nombre fini de composantes basiques distinctes 
(structures non etiquetees) qui permettent de proceder a une enumeration bijective (une 
approche tres rapidement impraticable, en particulier a cause du calcul des nombres 
d'automorphismes) . 

Remarque 2.2.13. Comme pour l'elagage, effondrer une hyperarete d'une chaine en un 
de ses sommets de degre > 2 ne change pas 1' exces de la structure. 

Definition 2.2.14. Une hyperarete speciale est une hyperarete contenant au moins 3 
sommets au moins de degre 2. 

Remarque 2.2.15. Une hyperarete speciale ne peut pas faire partie d'une chaine. 

Definition 2.2.16. Un sommet special est soit un sommet d'une hyperarete speciale soit 
un sommet de degre au moins 3 . 
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Ici, nous ne distinguerons dans les structures lisses que deux types de chaines : 

Definition 2.2.17. Une a-chaine est une chaine qui boucle sur un meme sommet. Une 
telle chaine admet au moins deux hyperaretes. 

Definition 2.2.18. Une (3-chaine est une chaine qui lie deux sommets speciaux. Une 
telle chaine admet au moins une hyper arete. 

Soit une composante B basique donnee avec les caracteristiques suivantes : 

r- d'exces i > 1 

ayant m sommets (non etiquetes) 
ayant s hyperaretes, necessairement s > [j^j + lj 
< - ayant c a a-chaines (2-19) 
ayant cp /3-chaines 
ayant p = c a + cp chaines 
de nombre (rationnel) d'automorphismes g . 

Alors les hypergraphes lisses (etiquetes) contenant la composante B a des effondrements 
pres de ses chaines, admettent la SGE 



1 F 

La lecture de cette equation est la suivante : 

1. Fixer un etiquetage des m sommets, c'est le facteur t m . 

2. Regrouper les sommets ainsi etiquetes dans la forme de la structure B, c'est le 
facteur 1 jg . 

3. Eventuellement, rallonger les p chaines dans l'ordre induite des etiquetages, c'est le 
facteur (1 - r(l))" p . 

Avec les caracteristiques de B , comme m = s(b — 1) — £, une reecriture de (|2.20p nous 
obtenons : 

Lemme 2.2.19. La SGE des composantes lisses (etiquetes) contenant la composante B 
basique d des effondrements pres de ses chaines : 

[(6-2)!]* r(t)« 

L'enumeration des composantes lisses complexes d'exces donne se reduit, ainsi grace 
a ce lemme, a la determination d'un ensemble de structures basiques permettant de 
les generer sans ambigiiite, c'est a dire qu'une composante ne peut etre generee qu'a 
partir d'une unique structure basique. Wright E.M obtient un ensemble de structures qui 
convient en considerant l'ensemble fini de toutes les structures basiques d'exces I n'ayant 
que des chaines de longueur minimale selon leur type, et precise ainsi la forme des SGEs 
des graphes connexes complexes lisses d'exces I . Ici, nous faisons le choix d'un autre 
ensemble de structures basiques pour montrer que 
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Theoreme 2.2.20. La SGE des composantes lisses d'exces £ est : 

avec les coefficients Ai p , des fractions rationnelles positives deb a coefficients entiers, 

r, = + 1J , (2.23) 

r®=(TW- (2 ' 24 » 

La forme (I2,22|) de la SGE des composantes complexes d'exces I donne une bonne 
intuition d'un terme principal qui contribue le plus au nombre asymptotique de ces struc- 
tures. En particulier, ce terme est lie a l'indice p = 2>i et, entre autres references, dans 
|19j , il est lie aux structures qualifiees de "clean" maximisant le nombre de chaines. 
Notons dans ce theoreme, le resultat combinatoire de la positivite des coefficients : 
la preuve est bijective et une preuve par induction sur (£, p) est restee hors de notre portee. 

D'abord une consequence immediate du lemme [2.2.191 est que 

Remarque 2.2.21. La SGE Hi des hypergraphes connexes lisses d'exces I s'ecrit sous 
la forme 

H £ (t) = i^, (2.25) 

avec 

0(t) = 1 - r(t) . (2.26) 

et fi un polynome de Laurent. 

Ensuite pour prouver le theoreme 12.2.201 . 

Lemme 2.2.22. Dans une composante complexe basique d'exces I, le nombre p de 
chaines est au plus 3£ soit p < 3£ . 

Preuve. II suffit de voir la preuve dans le cas d'une composante n'ayant que des chaines de 
longueur minimale. Soit done B une telle composante complexe avec les caracteristiques 
suivantes : 

r- d'exces £ > 1 

ayant so hyperaretes speciales 
ayant mo sommets speciaux 
ayant c a a-chaines (toutes de longueur 2) 
ayant eg /3-chaines (toutes de longueur une) 
ayant p = c a + cp chaines. 

Par definition de l'exces, 

m + 2(6 - l)c Q - c a + (b - 2) Cf 3 + £ = (b - l)(s + 2c a + cp) (2.28) 
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soit 

m -so(b-l)+£ = c a + cp = p. (2.29) 

Notons Bq l'hypergraphe induit par les sommets speciaux et soit £q son exces alors 
l'equation precedente s'ecrit : 

-£ + £=p. (2.30) 

La valeur de l'exces £q est minimale dans le cas oil Bq est une foret maximisant le nombre 
de composantes, Bq devant done se composer de sommets isoles de degre 3 ou d'hy- 
peraretes ayant exactement 3 sommets de degre 2 dans l'hypergraphe basique dont ils 
sont issus. Comme il existe un hypergraphe basique d'exces £ (voir la figure 12.41 ) ayant 
2£ sommets speciaux et sans aucune hyperarete speciale pour laquelle la valeur £q = — 2£ 
est minimale, nous deduisons 

-£ + £>2£ + £ = 3£>p. (2.31) 
Et cette borne superieure est atteinte. 




Fig. 2.4 - Structure d'un hypergraphe sans hyperarete speciale ayant le nombre 31 maxi- 
mal de chaines dans le cas d'exces £ . Seuls les sommets de degre > 3 y sont representes. 

Ainsi, dans (I2.22p , l'indice p est liee au nombre de chaines d'une structure basique. Quelle 
serait une interpretation du degre maximal du polynome fi de Laurent dans la remarque 
12.2.211 ? Dans l'equation (|2.22p , une interpretation intuitive du degre de r(t) est qu'il 
compte le nombre minimal d'hyperaretes dans une composante lisse. 

Lemme 2.2.23. Le degre maximum du polynome ft de Laurent, dans la remarque \2.2.21\ . 
est au plus 

Preuve. La demonstration passe par la decomposition des composantes lisses d'exces £ sui- 
vante : A une composante d'exces £ donne, supprimer la plus petite hyperarete dans l'ordre 
alphabetique induite par l'etiquetage, et obtenir ainsi des composantes eventuellement 
non lisses avec les caracteristiques (£i, k{) e'est a dire d'exces 1% et ayant hi sommets ap- 
partenant a la plus petite hyperarete supprimee. Des nombres de cycles nous deduisons 



£ + \ = Y,{^ + h) 

i 



(2.33) 
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Notons alors les divisions euclidiennes par (6 — 1) suivantes : 

e i + k i = q i (b-l)+r i (2.34) 

et 

£ + l = q(b-l)+r. (2.35) 

Ainsi, 

q(b-l)+r = Y^{Qi{b-l)+n} (2.36) 

i 

et _ 

q + l>^qi + l. (2.37) 

i 

En remarquant que le degre maximum de fe s'interprete comme le nombre minimum d'hy- 
peraretes, nous avons la demonstration par induction car le second membre de l'equation 
precedente maximise le degre de fi en prenant une composante ayant un nombre minimum 
d' hyper aretes . 

La forme du SGE Hi (I2.22|) suggere l'existence d'un ensemble de composantes basiques, 
d'exces i ayant s = ri, . . . , hyperaretes et dans laquelle seules s — rg chaines peuvent 

etre inserees, permettant alors de generer toutes les structures lisses. Le choix de prendre 
les composantes basiques n'ayant que des chaines de longueur minimale, bien que ce soit 
la justification immediate de la remarque 12.2.211 . ne permet pas d'aboutir au resultat car 
le nombre d'hyperaretes peut exceder strictement le nombre des chaines plus rg . 

La condition qu'une /3-chaine doit au moins contenir une hyperarete est motivee 
par la lecture plus facile des emplacements ou les chaines seront introduites. Aussi, en 
considerant Pensemble des composantes basiques n'ayant que des a-chaines qui elles 
memes sont de longueur 2 , toutes les composantes lisses peuvent etre generees sans 
ambiguite. Toutefois, la condition d'avoir autant d'hyperaretes que de chaines plus rg 
n'est pas une conditionalite respectee a priori dans la fagon canonique dont est faite la 
generation avec cet ensemble. Un plus apporte par la consideration d'un tel ensemble 
est qu'il garantit maintenant un nombre sumsamment peu eleve d'hyperaretes dans les 
composantes considerees d'apres les 2 lemmes qui suivent. 

Lemme 2.2.24. Pour une composante basique d'exces I ayant sq hyperaretes speciales 
et c a a-chaines 

2c a + s <3£ + l. (2.38) 
II existe une composante basique telle que I'egalite est atteinte. 

Preuve. Pour maximiser la quantite 2c a +so dans une composante basique d'exces £ donne, 
c a est a maximiser en priorite puis seulement apres sq . La relation qui lie le nombre des 
a-chaines et Pexces est c a < i + 1 et parmi les composantes basiques telles que c a = i + 1 
au plus so = t — 1 , nous deduisons ainsi le lemme. 

Une consequence directe de ce lemme est : 
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Lemme 2.2.25. Dans une composante basique d'exces £ n'ayant que des a-chaines toutes 
de longueur 2 , le nombre d'hyperaretes est au plus {*&£ + 1) . 

Preuve. (Theoreme I2.2.20D . Les composantes lisses peuvent etre generees sans ambigiiite 
par des composantes basiques ayant s = rg, . . . , 3£ + rg hyperaretes, c a a-chaines toutes 
de longueur 2 et cp /3-chaines toutes de longueur 1 et s — c a — cp marquages indiquant 
les emplacements ou d'autres chaines peuvent etre inserees. Notons que la quantite (s — 
c a — cp) doit etre positive et que les emplacements correspondent a des sommets ou une 
/3-chaine aurait pu etre effondree. Notons aussi que ces composantes basiques ne doivent 
pas etre toutes prises car dans ce cas, une meme composante lisse etiquetee peut etre 
obtenue par differentes composantes basiques. En bref, l'idee cle de la preuve : D'un cote, 
il y a l'ensemble de composantes basiques avec les chaines de longueur minimale pour leur 
type et, de l'autre cote, il y a l'ensemble de composantes basiques avec que des a-chaines 
de longueur 2 qui permettent de generer de maniere canonique toutes les composantes 
lisses. Entre les deux, il y a plusieurs choix d'ensemble de composantes basiques pouvant 
servir a generer toutes les composantes lisses . En particulier, il est possible de generer les 
composantes lisses a partir de composantes basiques ayant s = r^, . . . ,3£ + r£ hyperaretes 
en faisant grandir s — rg chaines aux emplacements des a-chaines toutes de longueur 2, 
des /3-chaines toutes de longueur 1 et des sommets sur lesquels une /3-chaine aurait pu 
etre effondree. 

Dans cette section, nous avons vu le nombre des hyperarbres enracines ainsi que 
celui des hypercycles. Quant au nombre des composantes complexes, nous pouvons le 
determiner du moins theoriquement car nous savons determiner les SGEs Hg o T{z) des 
composantes d'exces t donne : nous connaissons T{z) et nous savons comment obtenir 
Hg{t) . Cependant, la maniere vue jusqu'ici pour determiner ces SGEs est tres rapidement 
rebutante et tres difficile a automatiser, particulierement a cause de la determination des 
composantes basiques dont pour chacune il faudrait determiner son nombre d'automor- 
phismes ce qui nous parait peu concevable et surtout parce qu'il faut aussi pouvoir ga- 
rantir avoir considere toutes les composantes basiques necessaires, ni une de plus ni une 
de moins. Nous tirerons profit de la forme des SGEs Hi apprise dans cette section pour 
justifier comment les determiner recursivement dans la section qui suit. 

2.3 Enumeration recursive 

Ayant percu la puissance d'expression de l'outil que sont les SGEs, dans la presente 
section une etape supplementaire est franchie pour l'enumeration des composantes en 
etablissant une recurrence dont la resolution permet le calcul automatisable des SGEs Hi 
des composantes d'exces £ . 

Une breve apercue du contenu de cette section est la suivante : 

- elle commence par une vision recursive de l'enumeration des hyperarbres enracines ; 

- un outil pour la determination des coefficients des SGEs, la formule d'inversion de 
Lagrange, est presente avec une interpretation combinatoire ; 

- puis nous faisons une introduction aux SGEs bivariees qui servent pour une lecture 
facile de 
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- la decomposition des composantes d'exces t qui mene a une recurrence sur les SGEs 
univariees 

- ensuite la determination des Hi ou la resolution de la recurrence est presentee ; 

- et enfin nous terminons par la presentation de deux identites combinatoires per- 
mettant d'obtenir les sous la forme, permettant une lecture bijective, annoncee 
dans la section precedente. 

2.3.1 Enumeration recursive des hyperarbres enracines 

Les hyperarbres enracines sont des structures cles pour les hypergraphes : les SGEs des 
composantes connexes d'exces t s'expriment en fonction de la SGE T(z) des hyperarbres 
enracines. Le dictionnaire, introduit dans la section precedente, permet de traduire une 
decomposition naturelle des hyperarbres enracines comme etant composes d'une racine et 
d'un ensemble de groupes de (b — 1) hyperarbres enracines. Cette decomposition est une 
definition recursive des hyperarbres enracines. Effectivement, un hyperarbre se definit 
a partir d'hyperarbres de taille strictement plus petite. Cette definition recursive des 
hyperarbres enracines se lit dans une definition implicite de SGE via le dictionnaire. 

Definition 2.3.1. La SGE T(z) des hyperarbres enracines : 



La lecture de cette definition est la suivante : un hyperarbre enracine est 

- une racine, c'est le facteur z 

- et un ensemble (c'est la fonction exp ) non ordonne de groupes de (b— 1) hyperarbres 




(2.39) 



enracines (c'est I'exposant 





(2.40) 




Fig. 2.5 - Un hyperarbre enracine decompose. 



18 



CHAPITRE 2. ENUMERATION EXACTE 



Par cette definition, il est possible de determiner recursivement les coefficients [z n ] T{z) 
qui sont les nombres des hyperarbres enracines a un facteur n! pres. Evidement, cette fapon 
de calculer est a proscrire dans ce cas puisque les coefficients sont connus de maniere ex- 
plicite : l'expression mathematique (|2.5() pouvant etre retrouvee par la formule d'inversion 
de Lagrange via cette definition implicite de T(z) . 



2.3.2 Interpretation combinatoire de la formule d'inversion de Lagrange 

La formule d'inversion de Lagrange est un outil puissant pour faire de remuneration de 
structures arborescentes. La preuve de la formule d'inversion de Lagrange dans pEH page 
167] est une preuve analytique : l'extraction des coefficients des SGEs y est exprimee 
avec la formule integrate de Cauchy. Cette vision analytique de la formule d'inversion 
de Lagrange sera notre point de depart pour obtenir des equivalents asymptotiques des 
coefficients dans le chapitre 3 . Ici, nous adoptons un regard combinatoire dans la preuve 
du theoreme : 

Theoreme 2.3.2. Notons la SGE G 



00 + k 



G W = I>4' (2-41) 

k=0 

alors par la definition de la SGE T(z), la formule d'inversion de Lagrange s'ecrit 

[z n ] G o T(z) = ~ [t- 1 ] exp (nj^^j j-G(t) . (2.42) 



Preuve. Notons 



et 



Nous obtenons alors 



*W = exp ( jb^Ty. ' ' (2J;]! 



00 f m 

$(if= Y, Mn)— (2.44) 

m=0 



[ t n-l] ^ m nj_ G{t) \ = (2 . 45) 

E ^ n ^kL. = {2M) 

k,m>0\k+m=n—l 

= E ^n- fc -i(n) gfc+1 1 , (2.47) 

fc>0 v ' 
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(n-l)![t"- 1 ]|d>(t)«^G(t)}= (2.48) 

(n-1)! 
■fe!(ra-fc-l)! 



E^-iCnW .,^,, 1 !' , = (2.49) 



fc>0 



A:>0 



n — 1 



^24> n ~k-i{n)gk+i[ , ) = (2.50) 



4>n-k-i(n)g k+ i— — I I. (2.51) 

fc>0 

Observons que 

( ) -^[(6-1)!]--! [-(6-1)]!' (2 ' 52) 

soit 

n s (n — — 1)! 

0n - fc - l(n) = [(b-l)!]« s ! ' (2 ' 53) 

avec n = n(s) = s(b — 1) + k + 1 . A partir de (12.530 , nous identifions le nombre de forets 
de (k + 1) hyperarbres enracines a n sommets du theoreme 12.2.31 : 

n \ k + 1 , . , 

L'interpretation combinatoire est alors evidente en multipliant (|2.42p par n\ . Nous re- 
trouvons alors nos marques de la section precedente, pour enumerer les structures de 
taille n : 

- les gk+i structures lisses a (k + 1) sommets sont croisees, dans l'ordre croissante 
(uniquement definie) des etiquettes, avec 

- les forets de (k + 1) hyperarbres enracines a n sommets au nombre de 

n \ k + 1 , , , . 

k + i)— *-«<">■ (2 - 55) 

o 

Exemple 2.3.3. Nous retrouvons le nombre des hyperarbres enracines en prenant le plus 
simple des fonctions G , dans la formule d'inversion de Lagrange, c'est a dire la fonction 
identite. En effet, 

nl[z n ]T(z) = („-i)![t»-i]exp(n^^ (2.56) 



n 



71 

<"-%rwa- (2 ' 58) 
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avec n = n(s) = s(b — 1) + 1 . 

Exemple 2.3.4. Le nombre des forets de (k + 1) hyperarbres enracines de SGE 
est 

»»[*"] = (n - i)![ ^ l] (^ exp ( n (^)) (2 - 59) 

~ k\ J ^[( 6 -i)!]i/ J (2 - 60) 



(n — 1)! n s 

k\ [(&- l)l] s s\ ' 

avec n = n(s) = s(b — 1) + k + 1 , soit en accord avec le theoreme 12.2.31 



(2.61) 



, . „.T(z) k+1 ( n \„ , (n-jfe-1)! , , 

Exemple 2.3.5. De Pexemple precedent, nous deduisons une formule explicite du nombre 
des forets d'hyperarbres enracines : 

n—l Sk 

n\ exp(T(z)) = (n- 1)! ^ p , (2.63) 

avec n = n(s^) = 5^(6 — 1) + k + 1 pour tout A; = 1, . . . , n — 1 . 

Ainsi, par rapport a la section precedente, nous disposons maintenant d'un outil puissant 
cachant toute une machinerie bijective permettant de traduire, a partir des SGEs des 
structures lisses, les structures eventuellement non lisses correspondant. Ceci nous motive 
encore plus a automatiser la determination des SGEs des structures lisses. Pour y parvenir, 
dans un souci de clarte, nous presentons brievement les SGEs bivariees. 



2.3.3 Introduction aux series generatrices exponentielles bivariees 

A une sequence de nombre (a n ,m)n,me]N est associee une serie generatrice exponentielle 
bivariee 

A(w,z)= V a^ m w m —^. (2.64) 

^-^ TV. 

n,m£W 

Dans les hypergraphes, il y a deux notions essentielles a savoir la notion de sommet et 
la notion d'hyperarete. Soulignons que dans cette these, les series sont exponentielles en 
la variable z relative aux sommets qui sont etiquetes et elles sont ordinaires en la variable 
w , dans le cas de SGEs bivariees, relative aux hyperaretes. Si jusqu'ici, nous avons choisi 
de travailler uniquement avec des SGEs univariees, c'est parce que seul le parametre de 
la taille de la structure nous est pertinent et l'utilisation de la notion d'exces simplifie 
souvent de beaucoup les expressions des SGEs : c'est souligne dans [IB] par la richesse 
des resultats enumeratifs qui y sont obtenus. 
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L'introduction des SGEs bivariees nous permettra, en particulier, de souligner l'effet 
de decompositions dans les structures etudiees. Le dictionnaire de la section precedente 
s'enrichit du marquage d'une hyperarete ; les operations de marquage se traduisent en 
terme de derivation partielle : 



Operation sur les structures 


SGE correspondant 


Marquage de k sommets de A 
Marquage d'une hyperarete de A 


w 7h^A(w,z) 



Exemple 2.3.6. Pour enumerer les graphes etiquetes ayant 5 sommets et exactement 3 
aretes formant une clique d'ordre 3 nous pouvons recourir a la SGE bivariee : 



w°z 2 



W 3 Z 3 

~3T 



5\ w 3 z 5 



5! 



(2.65) 



Dans le membre gauche de cette equation : 

- le premier facteur correspond a la SGE des graphes ayant 2 sommets et sans arete, 

- le second facteur correspond a la SGE des graphes ayant 3 sommets et 3 aretes et 
le produit s'interprete a partir du dictionnaire pour justifier l'enumeration. 

La SGE bivariee de ces graphes avec une arete distinguee est 



5\ w 6 z 



w 



dw\2j 5! 



3^5 



5\ W 3 Z 5 

2 



Definition 2.3.7. La SGE bivariee des composantes d'exces I s'ecrit 



H e (w,z) 



w 



t I^H t oT[w 1 /^z) 



n w 



(2.66) 



(2.67) 



avec n = s(b — 1) 
et n sommets. 



et hg(s,ri) le nombre de composantes d'exces t ayant s hyperaretes 



Exemple 2.3.8. La SGE bivariee des hyperarbres enracines est 

T(wV(6-i)*) 

^1/(6-1) • 

Rappelons la definition implicite de la SGE T(z) 

'T{z) h - 1 



T{z) 



z exp 



(6-1)! 



Ainsi, 



dT(z) 



dz + Zj^-^dT(z)jexp 
m/ .dz T(z) b - 1 lm . . 

= T W7 + (Ri dr W- 



T(z 



,6-1 



(6-1)! 



(2.68) 

(2.69) 

(2.70) 
(2.71) 
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dT ( z ) = 1 dz 



(6-2)! 

Soit 

qui se lit comme l'equivalence des structures d'hyperarbre enracine ayant un second mar- 
quage (la racine peut etre eventuellement re-marquee) et une chaine orientee. Cette rela- 
tion (|2.73|) donne 

w dw { z) ~b-if, T( W m b -D z)b -i \ 

V (b-2)\ ) 

z dz { ] ~ fi nwvv>-i)z)<>-i \ ■ 

\ ( 6 - 2 ) ! / 

Nous en deduisons 

w-^H e (w,z) = (2.76) 
d 

b-i) 

H e oT(w 1 ^ b - 1 h)+ (2.78) 



W JL {wmb-D Hl T {w x '^z)) (2.77) 

V 1 (6-2)! / 

^oT(^'>z)| (2.80) 

^^or^ 1 )^ (2.82) 



J_*3 (wV^Ht o iW/M^ (2.83) 
— 1 oz V / 



soit la transcription d'un marquage d'une hyperarete, en terme de marquage d'un sommet 
pour les composantes d'exces t , suivante : 

w^Hi{w, z) = ^-j UH t (w, z) + ^H t {w, z)J . (2.84) 

La reside l'interet de classer les structures selon leur exces car les decompositions ex- 
primees en SGEs bivariees se ramenent alors a des SGEs univariees en reexprimant les 
SGEs avec un marquage d'hyperarete et en fixant w = 1 . 
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Exemple 2.3.9. Nous avons pour les hyperarbres enracines 

T(w l ^ b - l h) 



T(z). 



w=l 



Exemple 2.3.10. Et pour les composantes d'exces 



H e o T(z) 



(2.85) 



(2.86) 



w=l 



Exemple 2.3.11. Les composantes d'exces I ayant une hyperarete marquee sont enume- 
rees par la SGE univariee 



w=l 



^ (lH t o T(z) + z^-H e o T{z) 



(2.87) 



2.3.4 Recurrence des SGEs He 

Dans cette sous-section, l'exces I est fixe et a partir d'une decomposition des compo- 
santes d'exces t est etablie une recurrence des SGEs He ■ Pour decrire la decomposition, 
les SGEs bivariees sont d'abord utilisees car elles permettent une lecture directe des 
operations, sur les structures, qui y sont suggerees. Illustrons l'idee principale de la 
decomposition, menant a la recurrence, a partir de l'exemple de la figure [231 . 




Fig. 2.6 - La decomposition d'une composante sur un exemple. 

Dans cette figure, la figure a gauche represente une composante 4-uniforme ayant 14 
sommets et 5 hyperaretes dont une d'elle est marquee a savoir l'hyperarete {2,4,8, 12} . 
La figure a droite represente un hypergraphe 4-uniforme ayant aussi 14 sommets mais 
une hyperarete de moins que la composante a gauche. L'hypergraphe a droite n'est pas 
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connexe : il y a 3 composantes telles que chacune admet au moins un sommet marque. 
Les deux figures sont equivalentes l'une l'autre et representent la meme structure, seule le 
regard qui y est porte change : d'un cote c'est une hyperarete qui est marquee, de l'autre 
ce sont des sommets de cette hyperarete qui sont marques. II est clair que marquer une 
hyperarete dans une composante revient a marquer au moins un sommet dans chacune 
d'un certain nombre de composantes et inversement, en prenant soin qu'il y ait autant 
de sommets marques que necessaire pour former une hyperarete et qu'ils n'appartiennent 
pas tous deja a une meme hyperarete. 

II nous faut cependant faire attention pour un exces £ donne que les composantes avec 
chacune au moins un sommet marque se recombinent bien en un hypergraphe 6-uniforme 
d' exces £ . Pour cela nous faisons le lemme suivant : 

Lemme 2.3.12. Considerons une famille de composantes d' exces ji et ayant k{ sommets 
marques, —l<ji<£etl<ki<b. De plus, considerons que cette famille n'est pas celle 
constitute de seulement une composante avec b sommets marques tous appartenant deja 
a une meme hyperarete. Alors si une hyperarete est cree a partir des sommets marques 
des differentes composantes, Vhypergraphe obtenu est connexe et b-uniforme et d 'exces 
£ > —1 si et seulement si 



Preuve. II est clair que l'hypergraphe obtenu est connexe et qu'il est 6-uniforme si et 
seulement si le nombre total des sommets marques dans la famille consideree est egale a 



Notons par N le nombre de composantes dans la famille consideree et par n le nombre 
total de sommets de ces composantes. Dans la i-ieme composante, notons n« le nombre de 
ses sommets et Sj le nombre de ses hyperaretes. Alors, l'exces de la composante obtenue 
a partir de la famille vaut £ si et seulement si 



' J2 ki = b 



< 




(2.88) 



b : 




(2.89) 



N 




(2.90) 



Cette derniere equation est equivalente aux suivantes : 



N 




(2.91) 



N N N 

J2si{b-1) + J2 k i-J2 ni = e+1 ( 2 - 92 ) 

i=l i=l i=l 
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N 

(2.93) 



j>2({8i(b-i)-TH} + i*\ =e+i 



N 

J2tii + k i)=t+l- (2.94) 
i=l 

La decomposition, suggeree par le marquage d'une hyperarete dans les composantes 
d'exces £ > — 1 et illustree par la figure l2~6l . se traduit en terme de recurrence de SGEs 
bivariees dans le theoreme suivant : 

Theoreme 2.3.13. Pour tout £ > —1 , la SGE bivariee Hi des composantes d'exces I 
satis fait la relation : 

w^-H e {w,z)= (2.95) 



dw 



w 



U»Cy C ^ exp [ £ £ Cyc^^|^>, *) | + (2-96) 



fc=ii=-i 



- w -^-H e _ b+1 ( w ,z), (2.97) 



ou ff,- = si j < —2 . 



Preuve. La relation dans ce theoreme traduit en terme de SGEs bivariees la decomposition. 
II faut y lire : 

- dans le premier membre, la SGE bivariee des composantes d'exces £ avec une hy- 
perarete marquee (pour suppression) (wJj^Hi(w, z)) 

- dans le second membre, la SGE bivariee des families de composantes ayant k som- 
mets marques et d'exces j , lesquelles families, par creation d'une hyperarete definie 
par les sommets marques des composantes, produisent une composante 6-uniforme 
d'exces £ . Le terme avec un signe negatif contraint les marquages des b sommets a 
ne pas appartenir a une hyperarete deja existante (a ne pas recreer l'hyperarete, les 
structures considerees ici n'admettant pas d'hyperarete multiple) . L'extraction du 
coefficient de U b garantit d'avoir marque b sommets pour definir une hyperarete et 
l'extraction du coefficient de Cyc^ +1 garantit de former une composante d'exces £ 
d'apres le lemme precedent. Dans le premier terme de ce second membre, le facteur 
uu vient confirmer qu'il s'agit bien de creer une hyperarete de plus. 



Ce theoreme illustre la puissance d'expression des SGEs pour decrire les bijections ou des 
decompositions entre des structures combinatoires. En effet, dans ce theoreme est etablie 
une relation de recurrence entre les composantes d'exces £ d'un cote et les composantes 
d'exces strictement plus petit d'un autre. Cependant, Putilisation des SGEs bivariees 
pour cette relation de recurrence complique les notations car pour une valeur de l'exces £ 
donnee, la SGE He(w, z) bivariee se deduit de la SGE Hg o T(z) et inversement bien sur. 
En effet, rappelons (12^861) et (p^7j) : 
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a partir de la SGE bivariee se deduit la SGE univariee (il suffit de fixe la valeur de 
la variable uu = 1) 



= H e o T{z) , (2.98) 

UI = 1 



- a partir de la SGE univariee se deduit la SGE bivariee pour un exces donne (par 
definition de l'exces, le nombre de sommets est lie a celui des hyperaretes) 

H e {w, z) = wW-^Hz o T{w x t^z) ■ (2.99) 

Nous pouvons transcrire le theoreme 12.3.131 en terme de SGE univariee, les derivees 
partielles deviennent des derivees droites : 

Theoreme 2.3.14. La SGE H^o T(z) des composantes d'exces i satisfait la relation 

' eH e oT(z)+z—H e oT(z))= (2.100) 



IV dz 

£/ 6 Cyc^] exp [ £ {I Q^^^i ° T(») I + (2-101) 
\k=ij=-l 

' ' 1 + l)i^_& + i o T(z) + z-^-Ht-w o T{z)) . (2.102) 



6-1 V dz' 

Preuve. Avant de fixer w = 1 , nous observons que le premier membre de l'equation du 
theoreme 12.3.131 se reecrit (rappel de (|2.84p ) 

w^H t {w, z) = ^ (eHt(v, *) + z-j^M™, «)J • ( 2 -l° 3 ) 

Et ainsi, nous deduisons le theoreme. 

Par ce theoreme, nous disposons d'une recurrence des SGEs HgoT(z) sous la forme d'une 
equation differentielle d'ordre un (dans la mesure ou les termes du second membre qui 
apparaissent avec Hj oil j < £ ont ete explicites dans une resolution anterieure) . 

Nous avons affirme sans justification, dans la remarque 12.2.91 . que la SGE des hyper- 
arbres non enracines s'ecrit sous la forme H-\ o T{z) , soit en fonction de la SGE T(z) 
des hyperarbres enracines. La justification est ici faite par l'application de ce theoreme 
pour la valeur I = — 1 , soit le cas de la SGE des hyperarbres non enracines : 

Theoreme 2.3.15. La SGE H-\ °T(z) des hyperarbres non enracines est telle que 

H- 1 (t)=t-t=^. (2.104) 

Cette ecriture est valable dans le cas ou b = 2 , c'est a dire dans le cas des arbres non 
enracines. 

L 'ecriture du theoreme 12.3. 14l se simp line en une equation differentielle en la "variable" 
T(z) : 
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Theoreme 2.3.16. Pour £ > 0, notons qu'il existe une fonction telle que 

J £ oT(z)= (2.105) 

,(Uz) k d k 



U b Cyc i+1 
1 



k=ij=-i 



M d^ oT(z) l + (2 - 106) 



b-i 



1 + l)Hi~ b+1 o T(z) + z— He-b+i o T(z] 



Ainsi, 



6-1 



£H e (t) + t—H e (t) ) =J e (t). 



(2.107) 



(2.108) 



Preuve. Dans l'equation ()2.100p , en regroupant les termes avec l'indice, relative a l'exces, 
egale a I dans le premier membre nous obtenons : 
- d'un cote 



1 



6-1 



soit par (gTTgp 



T(z) b ~ 1 \ d 
£H e o T(z) + ( 1 - } ' ^. ) z-^H e o T{z. 



1 



2)! J dz' 



(2.109) 



de l'autre, l'intervalle de l'indice j devient j = — 1, 
l'exp , 



-(£H e oT(z)+T(z)H/ oT(z)) (2.110) 

, £ — 1 dans l'exposant de 



U b Cyc £+1 



exp 



fc=ii=-i 



6-1 



i + i)i^_ 6+ i o t(z) + z— o r(/ 



(2.111) 



(2.112) 



qui vaut (|2.110|) et comme ce dernier s'exprime en T{z) , nous deduisons l'existence 
de Ji telle que 



Jt ° T(z) 



(2.113) 



U b Cyc i+1 
1 



cxp 



.fc=lj=-l 



6-1 



(£ - 1 + l)if*_ 6+ i o T(z) + z—H^ b+1 o T(z) ) . (2.115) 



En substituant T{z) par la variable t , on obtient l'equation differentielle du theoreme : 



1 



6-1 



£H e (t)+t—H e (t) ) =J e (t). 



(2.116) 

❖ 
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Remarque 2.3.17. Une variante (plus bijective) de la preuve de ce theoreme est que 
la decomposition decrite dans Ji se comprend grace aux structures lisses : les structures 
obtenues par le marquage de sommets dans une composante peuvent etre classees en des 
composantes lisses. Bref, z k -^^Hj o T(z) s'exprime en fonction de T(z) . 



2.3.5 Resolution de la recurrence des SGEs Hi 

Ainsi, nous disposons, dans le theoreme 12.3.161 . d'une recurrence des SGEs Hi des 
composantes lisses d'exces I sous la forme d'une equation differentielle d'ordre un sui- 
vante : 

(tH t (t) + t^Ht(t) \ = Je(t) , (2.117) 

avec Jg defini a partir de (|2. 105j) . Dans cette sous-section, nous resolvons ces equations 
differentielles iterativement en l'exces I . A partir de la peut se faire l'automatisation du 
calcul des SGEs Hi des composantes lisses d'exces t > . La solution est uniquement 
determinee en considerant la condition initiale 



H e (t)\ t=0 = 0. (2.118) 
Adoptons la notation, rencontree dans la remarque 12.2.211 . 



$(t) = 1 - r(t) . (2.119) 

Remarque 2.3.18. Le choix de fixer cette notation, bien que 6 n'admette pas de lecture 
combinatoire directe, s'est impose a partir de Particle [22] puis, a posteriori, parce qu'il est 
commode de manipuler les polynomes de Laurent f# de la remarque 12. 2. 21 1 pour definir les 
SGEs Hi . Notons tout de meme que les series, en t , l/0(t) et (1 — (t)) admettent chacune 
une lecture combinatoire a partir desquelles les SGEs exprimees avec 0(t) peuvent etre 
comprises. 

Rappelons la definition des fy. 

Definition 2.3.19. fa est un polynome de Laurent defini par l'intermediaire du SGE 
H e : 

Hi(t) = l^. (2.120) 

Remarque 2.3.20. L'equation homogene associee a (|2.117p admet pour solution Cste/i^ . 
Ce qui nous suggere de trouver une solution sous la forme 



H e (t) = 




(2.121) 



Cette forme decoule aussi de ([2. 120j) . 



Par le lemmeEEE], 
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Remarque 2.3.21. Les composantes d'exces j ayant k sommets marques admettent une 
SGE qui s'ecrit sous la forme 



avec fjf. un polynome de Laurent. En particulier, fj$ = fj . 



(2.122) 



Exemple 2.3.22. Pour les hyperarbres, la forme de l'ecriture de la definition 12.3.191 est 
valide : 

'b-l + 9(ty 



soit 



ff_i(t) = t 



f-i o 0(t) 



b-l + 0(t) 



(2.123) 
(2.124) 



C'est le seul des fj de degre minimum positif (ici, ce degre est nul) 



Convenons de noter les polynomes de Laurent fjk avec la variable x , alors en evaluant 
^{-^kHj o T(z)} par derivation de fonctions composees sachant (|2.73|) 



hk+i(x) = ~(b - 1)^ + (6 - l)f jk '(x) - j 



fjk(x) 



X 



kfjk(x) 



Remarque 2.3.23. La determination recursive de fe ■ Comme 
Jt(t) = 



b £-1 



U°Cyc 



exp 



h=ij=-i 



V 



+ 



t e-b+i 



— f(p-i , ^ /<-hi,o°^) , fe-b+i,i Q 0(t) 
-iy L ~*~ L ) ft-b+i + t e ~ b+1 

U b Cyc £+1 



(2.125) 

(2.126) 
(2.127) 

(2.128) 



b l-l 

ex P ( £ £ Cyci +k U k f jk o 0(t) ] + (2.129) 
k=ij=-i 



1 ((e-l + l)f e - b+1 , o o6(t) + f e _ b+1A o0(t)) 
~ { b-V tf- b + l 



(2.130) 



Et comme en notant Hi sous la forme (|2.12ip , l'equation differentielle (|2.117p devient 

1 9&{t) 



soit 



b-l—- Je{th 
g/(t) = (b-l)t e ~ 1 J e (t). 



(2.131) 
(2.132) 
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Nous aurons 
9t(t) 



(2.133) 

(b e-i \ 
£ £ ^ +k U k f 3 k o 9{t) + (2.134) 
fc=ij=-i / 
i6 _ 2 ((*-! + l)//-H-i,o o 0(t) + //-6+1.1 o 



-{b-l)t b 



= -(6-2)! 



-(6-2)! 



dt 



U h Cyc i+1 



6-1 < 2 - 1SB > 
exp Cyc? +k U k f jk oO(t) ] + (2.136) 



d 0(t) + l)ft- b+ i, Q 0(t) + / f - 6+ i,i o gft)) 



dt 



car 



d0(t) 



6-1 
(6 - l)t b - 2 



dt. 



(6-2)! 

En identifiant 9{t) = 9 , 
g/(t)dt = 

(b e-i 
E E Vy<? +k U k fjm \d9+ (2.140) 
fc=lj=-l 



(2.137) 
(2.138) 

(2.139) 



. (6 _ 2) , ((^ - 1 + i)Wg) + wg» d 6 , 



6-1 



(2.141) 



Ce qui, par integration, en prenant la solution qui verifie fe(l) = , determine 

ft ° 0(t) = g e (t) . (2.142) 

Nous pouvons maintenant resumer la resolution de la recurrence dans le theoreme 
suivant : 

Theoreme 2.3.24. La SGE Hi des composantes lisses d'exces t > est determinee 
recursivement sous la forme 



avec 



ou 



ft,o( x ) = M x ) 



J<(«)d«, 



J e (u) = 



2)! [u b Cyc e+1 ] exp ( E E Cyc? +k U k f jk (u) ] + 
\ k =ij=-i J 

,{(£-! + l)WoH + fe-b+i,i(u)) 



-(6-2)!- 



6-1 



(2.143) 
(2.144) 

(2.145) 
(2.146) 

(2.147) 
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avec 

f hk+l (x) = -(b - + ( b - l)f jk '(x) - j^^- - kf jk (x) . (2.148) 

Notons en particulier que /_i i i(x) = 1 . 

Par ce theoreme, il est naturel de determiner les fi sous forme de polynome de Laurent, 
une forme simple qui n'admet pas de lecture combinatoire a premiere vue. Dans la sous- 
section suivante, nous montrons comment automatiser la determination de Hi sous la 
forme (|2.22p avec laquelle nous disposons d'une lecture combinatoire. 



2.3.6 Mise en forme des SGEs H t 

Avant de donner les premiers exemples des SGEs Hi , nous donnons deux identites 
combinatoires qui permettent de mettre les SGEs Hi sous la forme de l'equation (|2.22p . 

Lemme 2.3.25. Pour j, a G IN* (done j + a > 0) , 

^E(T)^+£(^-T>->"-*- 

i=0 x ' i=0 v J 7 

Lemme 2.3.26. Si a — j > alors 

03 = E ( a "_ J "_T x ) c 1 - e ) a ~ l ■ ( 2 - 15 °) 

oil si k £ IN et t £ 7L alors 

't\ ( t \ t(t — l)---(t — k + l) 



k \t-kj k\ 



(2.151) 



Les preuves de ces deux lemmes sont donnees en annexe. 

Par ces deux identites presentees dans ces deux lemmes, les deux formes des SGEs Hg 
se deduisent l'une de l'autre. Rappelons ici ces deux formes : 

• avec fi un polynome de Laurent, nous avons une forme pratique, 

H t (t) = (2.152) 

• a partir de (|2.22p , nous avons une forme combinatoire, 

(1- , fl-0(t"* 

p=0 

sachant r{t) = 1 - 0(t) = t b ~ l /{b - 2)! . 



ffi ( t)= ™£I E ^(i^p) (2 . 153 ) 
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La forme pratique se deduit de la forme combinatoire directement en developpant l'ex- 
pression de la forme combinatoire : 

= £ E CiWWW, (2-154) 

j=-3i 

avec 

r * = + !J • (2- 155 ) 

Et la forme combinatoire se deduit de la forme pratique en utilisant les lemmes avec la 
valeur a = ri pour reexprimer les 9{t) 3 . 

Remarque 2.3.27. Dans ce qui suit, dans le cas des graphes, c'est a dire si b = 2 , nous 
adoptons la notation Wi = Hg de Wright pour les series generatrices des composantes 
d'exces i . 

Definition 2.3.28. Notons X la SGE des chaines non vide 

m- 1 -^- (2.156) 
Nous sommes maintenant en mesure de donner les SGEs Hi des composantes d'exces 

t : 

Theoreme 2.3.29. La SGE H-\ o T(z) des hyperarbres non enracines est telle que : 

H- lW =t-^=*(*^)). (-57) 

Theoreme 2.3.30. La SGE Hq o T(z) des hypercycles est telle que : 
- sib = 2 (pour les graphes) 

l-6(t) (I -Of 



Wo{t) = -ln(^))-i^-ii^, (2.158) 



t t 2 



W (t) = -In(vT^t)--- -. (2.159) 

- si b > 3 

flo(t) = " In (VW)) ~ — |^ • (2-160) 

Theoreme 2.3.31. La SGE Hi oT(z) c/es composantes d'exces 1 es£ £e//e que : 

- sib = 2 (pour les graphes) 



si b = 3 



(l-9(t)) 2 f5 v , .3 , 19 w+ . 2 



= " r^- ( g*(t) d + + QX(t) ) . (2.162) 
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- si 6 > 4 

\-e{t) ( h{b-ifx(tf , 

Hl{t) = ^^{ 24 + 

(76-12)(6-l)X(t) 2 (b-2) 2 X(t) 
24 12 

Theoreme 2.3.32. SGE H2 o T(z) des composantes d'exces 2 es£ £e//e <j"ue : 

- sib = 2 (pour les graphes) 

5t 8 

55t 7 73t 6 3i 5 t A 

+ ^7777 7TT + "77777 777 + T77 777 + 



si 6 = 3 



48(1 -t) 5 48(1 -t) 4 4(1 -i) 3 24(1 - t) 2 



sib = 4 



t 2 V 3 



+ ^X(t) 4 + 118X(t) s + ^X(i) 2 + H X (f) 



si 6 > 5 



g ,( t ) = iz|W(w-'W + 

w t 2 V 16 

5(116 -17)(6-l) 3 X(t) 5 
48 

(26- 3)(386 - 65)(6 -l) 2 X(t) 4 
+ 48 + 

(6 - 1)(486 3 - 2446 2 + 4116 - 229)X(f) 3 
+ 48 + 

(6 - 2) 2 (136 2 - 446 + 35)X(t) 2 (6 - 2) 2 (6 - 3) 2 X(t) 
+ 48 + 48 

Theoreme 2.3.33. La SGE H3 o T(z) des composantes d'exces 3 est telle que : 
- sib = 2 (pour les graphes) 

1105t n 395t 10 1513lg 

aW ~ 1152(1 - t) 9 + 72(1 -tf + 1152(1 -t) 7 + 

2399t 8 8303i 7 557t 6 3t 5 

+ TTT77 7777 + 7777777 7777 + 77777 777 + 



144(1 -t) 6 720(1 -t) 5 144(1 -t) 4 8(1 -if 
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si b = 3 



si b 



si b = 5 



si b > 6 



m) ^l^)l^ x{tf + ^ x{tf+ (2 , 78) 
+ ™p x(()7+ 3|p x(<)6 + 1^21 x(()5+ (2179) 

7397 w . , 4 42701 v , ,» 4553 w , 2 41 \ 



H, (t) = llz|W)!(^ (t) . + ^ (t) . + (2 . 18 1) 
217431 , + 713511 x(()6 + 14467311 + 



16 w 32 w 640 

+ ^X W + ^X, ( )3 + ^X, t)2+ (2.83) 



m) = ilzMl! + ™»° * (t) . + (2.185) 

+ 8272K( t )« + ^-Xiff 4- (2.186) 
12565331 x(()4 + 4851107 + 46821 x(()2 + 



360 v 7 480 w 32 

1 

4~ 



2Q1 

+— X(t) ) . (2.188) 



g , W = l^( "«(>-W)' + ( , 189) 

5(12596 - 1922)(6 -l) 5 X(r.) 8 
+ 1152 + l2 ' iyUJ 

(6 - 1) 4 (151066 2 - 471086 + 36643)X(t) 7 

+ hm + (2 - 191) 

(6 - 1) 3 (98676 3 - 473686 2 + 755 926 - 40080)X(rQ 6 



576 

(6 - 1) 2 (755 2 96 4 - 49 9 5 6 46 3 + 123 4 8 6 06 2 - 13511526 + 551736)X(f.) 5 

5760 + 
(6 - 1)(337916 5 - 2917426 4 + 10029916 3 - 17150006 2 + 14570886 - 491520)X(t) 4 

5760 



(2.193) 
(2.194) 

(b - 2) 2 (20566 4 - 144246 3 + 373536 2 - 421486 + 17404)X 3 

H h (2.195) 

1440 V ; 

(6 - 2) 2 (566 4 - 5146 3 + 17536 2 - 26236 + 1444)X(t) 2 

+ 360 + (2 ' 196) 

(6-2) 2 (6-3) 2 (6-4) 2 X(f.) \ ^ (2 m) 
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Remarque 2.3.34. Les SGEs Hi des composantes lisses d'exces I pour b > £ + 3 ad- 
mettent une meme expression fonction de b . 



Lemme 2.3.35. Les fjk de la remarque \2.3.21\ sont des polynomes de Laurent de degre 
minimum egal a 

-(3j + 2fc)-l si (j,k) = (—1, 1) 
-(3j + 2fc) /(-!,!)• 



(2.198) 



Preuve. La preuve est immediate par la recurrence des du theoreme 12.3.241 . 







Lemme 2.3.36. Considerons les families de composantes d'exces ji et ayant ki sommets 
marques, — 1 < ji < £ — 1 et 1 < ki < b, composantes qui se recombinent en composante 
b-uniforme d'exces £ et sans hyperarete multiple. Seules les families ayant exactement 
(b — 2) hyperarbres enracines suffisent a la determination du coefficient du terme de degre 
minimum —3£ de fi . Et ces families avec celles ayant exactement (b — 3) hyperarbres 
enracines suffisent a la determination du coefficient du terme de degre (—3£ + 1) de fi . 

Preuve. Notons par m_i 5 i le nombre de hyperarbres d'une famille fixee. Alors le degre 
minimum des termes (dans ft) liees a la famille consideree est 



- ^(3ji + 2fcj) - m_ M - l = _(3£ + 3) + (6 - m-i,i) 

i 



"I, 



(2.199) 



d'apres le lemme precedent et la decomposition relatee dans l'expression j du theoreme 
12X241 . 

Ce dernier lemme nous permet d'expliciter les deux premiers coefficients de fi a savoir 
celui du terme de degre respectivement — 3£ puis (—3^+1) pour les composantes complexes 
{£ > 1) . Notre motivation pour mieux inspecter ces coefficients vient de leur importance 
dans la litterature pour faire de remuneration asymptotique comme souligne dans |28l 
page 42] : le premier coefficient est lie a ce qui est aussi denomme coefficient de Wright- 
Stepanove (voir ou encore coefficient de Wright-Takacs-Louchard et son expression 
asymptotique a fait l'objet d'etude [6l[7] permettant alors d'obtenir des enumerations de 
composantes d'exces tendant vers l'infini avec le nombre de sommets [24j . 



(2.200) 



Theoreme 2.3.37. La SGE des composantes d'exces £ > 1 est telle que : 



H e (t) 



1 X e (b - l) 2e _ ( K< -^(b- 2 ))(ft-l) 
t e 3£6(t)^ {3£ - l)^(t)^- 1 



2£-l 



j=-3l+2 



avec 



1 

2 



1 1 ^ l 

-Vi(3£-i)+-x;^- 1 - 



(2.201) 

pour £ = 1,2... (2.202) 



p=0 
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vi = h + h Xo (2 - 203) 



soit fio tel que : 



+ *s\ P \e-2-s-p (2.204) 

s=0 p=0 

£—2 

+ 2 A *-i + ^(3p + 2)A P A^ 2 _ P (2.205) 

+-(3£ - 4)(3£ - 2)A^_ 2 powr £ = 2, 3, . . . , (2.206) 
6 

Ho = b-l (2.207) 



alors 



1 

= -{{3£-2)M_ 1 + (3b£-b-2£)\ e _ 1 )+Y,» P ^-i-p, (2-208) 



2 

era /i£ est a calculer au besoin 



p=0 



2 

Ht = m -Mb -2) + Xe{b- -) pour* =1,2... (2.209) 

(Notons que Kg = ne(b)) . 

Dans ce chapitre, nous avons obtenu des enumerations explicites des hyperarbres 
enracines (done aussi de ceux non enracines) , des forets de (k + 1) hyperarbres enracines 
(done par deduction des forets d'hyperarbres enracines) et des hypercycles. La complexite 
rencontree dans les expressions des enumerations explicites nous a mene a expliciter non 
pas le nombre meme des composantes d'exces donne, mais leur serie generatrice via la 
notion des composantes lisses, et permettre ainsi de proceder a remuneration en utilisant 
la formule d'inversion de Lagrange. Si la determination des SGEs a pu etre automatisee 
grace entre autre a la grande expressivite des SGEs pour decrire des decompositions 
des structures a enumerer, la complexite des expressions explicites nous force a preferer 
des equivalents asymptotiques. L 'enumeration asymptotique est traitee dans le chapitre 
suivant. 



Chapitre 3 



Enumeration asymptotique 



Dans ce chapitre, notre but est de fournir un equivalent asymptotique du nombre de 
composantes d'exces t quand la taille, c'est a dire le nombre de sommets, est grande. 
Notre motivation vient, d'une part, de la difficulte du calcul du nombre exact pour des 
composantes de grande taille, et d'autre part, de la connaissance des SGEs des compo- 
santes lisses permettant de faire de l'estimation asymptotique via la formule d'inversion 
de Lagrange, la formule integrale de Cauchy en analyse complexe et, ici, notamment par 
la methode du point col pour obtenir nos resultats. Les notions relatives a l'utilisation des 
series generatrices pour remuneration asymptotique apparaissent dans [26J , en particulier 
ici les resultats asymptotiques sont obtenus via la methode du point col - voir [151 E! pour 
plus de precisions sur les methodes asymptotiques. 

Le plan de ce chapitre est le suivant : nous commengons par l'enumeration asympto- 
tique des hyperarbres enracines, puis enchainons avec celle des hypercycles et demontrons 
un encadrement des coefficients [z n ] Hi o T{z) , t > 1 , pour finalement apres etablir le 
nombre asymptotique des composantes complexes d'exces donne. 



Dans cette section, etant vue l'importance des hyperarbres enracines pour l'enumer- 
ation exacte des hypergraphes, nous procedons a l'enumeration asymptotique de ces com- 
posantes les plus simples. La methode du point col, pour faire de l'enumeration asympto- 
tique, est ici illustree et validee parce que le resultat se verffie aussi par une autre preuve : 
via la formule de Stirling, connaissant l'expression du nombre des hyperarbres enracines 
a un nombre de sommets donne. 

Nous commengons par un exemple, en fixant la valeur de b = 3 , et determinons 
le developpement asymptotique des coefficients de la SGE des hyperarbres enracines 3- 
uniformes. 

Proposition 3.1.1. Le nombre des hyperarbres 3-uniformes enracines ayant (2s + 1) 
sommets est determine a partir du coefficient de z 2s+1 du developpement de la SGE 



3.1 Enumeration asymptotique des hyperarbres 
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Ce coefficient admet le developpement asymptotique, quand (s — ► oo) , suivant 

I 2 1 r «« " WS5 " 24^ + iW " 414720? + ' ' J ' (3 ' 2> 

Preuve.(Vi& Stirling) Par le corollaire 12.2,41 . le nombre des hyperarbres 3-uniformes en- 
racines ayant s hyperaretes et n = n(s) = 2s + 1 sommets est 

n\ s i n! . . . i , 

-,n s = r(2s + l) s_1 . (3.3) 

2 s s! 2 s sr ' y ' 

Le coefficient de z 2s+1 de la SGE T^(z) est done 

l* 2S+1 ] W = • (3-4) 



Par l'expansion de Stirling nous obtenons 

2sV2~KS (! + T2i + 288^ ~ 51840s 3 ~~ ' " 



fl _|_ J_\s-l p s I 



Comme 

(i+ 2s j " 6 V 8s + 384s 2 1024 S 3+---J W 
nous obtenons l'expansion asymptotique du coefficient z 2s+1 de 13(2;) suivant 

r 2,+n nn , x e s+1 / 2 / 17 481 96133 \ , „ s 

I 2 1 " W55 I' - 547 + 1152? - 4T4720? + ' ' j ' < 3 ' 7 > 

<> 

Nous retrouvons une version de cette proposition 13.1.11 via la formule d'inversion de 
Lagrange et via la formule integrale de Cauchy : 

Proposition 3.1.2. Le nombre des hyperarbres 3-uniformes enracines ayant (2s + 1) 
sommets est determine a partir du coefficient de z 2s+1 du developpement de la SGE T%(z) 
definie par V equation (|3.1|) . Ce coefficient est, pour (s — ► 00) , tel que 

Preiive. (Methode du point col) Par la definition implicite p. II) de la SGE 13(21) , la formule 
d'inversion de Lagrange donne : 

[z 2s+1 ] T 3 (z) = -L^ [t 2s ] exp ((s + l/2)t 2 ) . (3.9) 

Dans le second membre de cette equation, en extrayant le coefficient de t 2s par la formule 
integrale de Cauchy, nous trouvons 
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ou le contour integral encercle l'origine du plan complexe. Ainsi, le coefficient de la SGE 
s'ecrit 

[^ 2S+1 ] Uz) = —L— ^gp ((« + l/2)(t 2 - ln(t 2 ))) dt. (3.11) 

Le parametre s , qui denombre les hyperaretes, est destine a tendre vers l'infini. Nous 
soupgonnons dans Pequation precedente un contexte dans lequel la methode du point col 
peut etre appliquee. Notons que la derivee de l'exposant dans Pintegrande de Pequation 
s'annule en deux points t = ±1 : 

±(a + l/2)(t 2 - ln(i 2 )) = 2(s + l/2)(t - - f ) . (3.12) 

Nous prenons alors comme contour integral le carre centre en , aux cotes paralleles 
aux axes et de diagonale 2y/2. Ce carre passe par les points cols t = ±1 . Dans la 
figure 13.11 . en representant quand t prend ses valeurs sur ce contour la valeur de la 
partie reelle de (t 2 — ln(t 2 )) , nous voyons que quand s est grand la contribution, dans 
le module de l'integrande, des piques aux deux points cols est tres accentuee, alors deux 
voisinages limites autour de chacun des points cols suffiront pour capturer l'asymptotique 
du coefficient recherche. Notons Pequation parametrique du carre 7 



1 








1, v, 1**2-v**2-ln(1**2+v**2) - 

-1, v, 1**2-v**2-ln(1**2+v**2) 

u, 1, u**2-1"2-ln(u**2+1**2) 
u, -1, u**2-1**2-ln(u**2+1**2) 

1,v, -1 

-1,v, -1 

u, 1,-1 - 

u, -1,-1 


0.5 





















-0.5 










-1 
-1 








1.5 








^5— — 




Fig. 3.1 


- La valeur de 5ft (i 2 - ln(t 2 )) 


sur le contour carre choisi. 








71 : t = 1 + iv , 


v]e [-1,1] 






7 : < 


7_i : t = -1 + iv , 
7i : u + i, 
_ 7_i : u + i , 


f i€ [-1,1] 
«i€ [-1,1] 
n|e [-1,1]. 



(3.13) 
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Soit 71 un voisinage de t = 1 du chemin 71 . Pour avoir l'ordre de grandeur, nous appro- 
chons, dans l'integrande, l'exposant par 

»((a + l/2)(t 2 - ln(t 2 ))) = K(( s + 1/2)(1 + 2(* - l) 2 + 0((* - l) 4 ))) . (3.14) 

Nous consider ons alors l'integrale 

h = 2vk(2s + 1) 1 6XP ((g + 1/2)(t ' ~ dt (3 ' 15) 

= l^+l) X 6XP ^ ((2S + 1){{t - 1)2 + - 1)4))) ) dt • (3 - 16) 



En effet, avec le changement de variable 

f t= 1 + 



y / 2i+I 



dt = .'. dr, 

V2s+1 



(3.17) 



nous avons 



exp (s + l/2) r exp / 2 + 0(7 ^ ))dr (3 . 18) 



27r(2a + 1)3/2 7_ K< ^ v (2s + 1) 

exp(s + l/2) [*' _ 2Wl , r4 

2tt(2s + I) 3 / 2 ' ' 
exp(s + 1/2) 



f[ exp(-r 2 )(l + 0(-^ T - T y))dr (3.19) 

W { H exp (_r2) dr + II exp (_r2) 0( (2^TT) )dr } (3 ' 20) 



2vr(2s + 
r etant une variable reelle, 

„4 



. ex P (- 2 ) (p JTly )d r = / W r (3,1) 



5 



0( (2jtl)»- <3 - 22) 



Ainsi, sur la portion 71 du chemin 71 , approcher l'exposant permet d'avoir 

5 



Comme 



f'QO f'QO 

/ exp(-r 2 )dr < / exp(-K s r)dr (3.24) 

1 r°° 

< — / exp(-r)dr (3.25) 
= Q( eXp( -^ 2) ), (3.26) 
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en prenant k s = (2s + l) 1 / 6 pour definir 71 , nous trouvons sur ce chemin 



T exp(s + 1/2) ( f°° . 2 , 1 
Nous y reconnaissons l'integrale de Gauss et pouvons ecrire que 



1)1/6- 



exp(s + 1/2) 
20F(2s + I) 3 / 2 



l + 0( 



1 



(2s + l)V6- 



(3.27) 



(3.28) 



Des lignes de preuve similaires s'appliquent pour la portion 7_i du contour. Sur cette 
portion, en prenant le voisinage 7_i centre sur le point col — 1 et de longueur deux fois 
n s , qui peut etre pris egale a (2s + l) 1 / 3 , l'exposant de l'integrande peut etre approche 
et nous obtenons 



/-1 = 



exp(( S + l/2)(t 2 -ln(t 2 )))dt 

... 5 

exp (-r 2 ) dr + 0(- 

i 

l + 0(- 



2m(2 

exp(s + 1/2) 
2vr(2s + I) 3 / 2 
exp(s + 1/2) 



•(2a + l) 



2^F(2s + I) 3 / 2 



1 



(3.29) 
(3.30) 
(3.31) 



(2s + l)V6' 
Soit la portion de 71 

j[ : t = 1 + iv, v]£ [k s /V2s + 1,1]. (3.32) 
Nous approximons l'exposant de l'integrande sur cette portion au point t = \+in s / 'y/2s + 1 . 
1 



2ivr(2s + 1) 



exp((s + l/2)(t 2 -ln(t 2 )))dt 



7i 



< 



1 



2vr(2s + 1) J Ks /^+i 
exp(s + l/2) f 1 



exp ((s + 1/2)(1 - v 2 - ln(l + t> 2 ))) d'. 



2tt(2s + 1) J Ks/y ^+T 
exp(s + l/2) rV2l+I 



[ exp((s + l/2)(-2v 2 + 0(v 4 )))dv = 



2tt(2s + I) 3 / 2 J K 



exp -r 2 + 0{- 



(2s + l) 

exp(s + 1/2) Z-^+i - 4 
2vr(2s + 

exp(s + 1/2) exp(-K s 2 ) 
2vr(2s + l) 3 / 2 1 



) dr = 



W-L ^(- r2 )< 1 + °<(2JTi)» d ' 



(3.33) 
(3.34) 
(3.35) 
(3.36) 
(3.37) 
(3.38) 



L'ordre de grandeur ci-dessus est valide pour k s = (2s + l) 1 / 6 et il souli gne que, sur la 
portion de 71 , la contribution de l'integrale est exponentiellement petite en dehors de 71 . 
Soit la portion de 7, 

j' i :t = u + i, u[e [0,1]. (3.39) 
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Nous approximons l'exposant de l'integrande sur cette portion au point t = 1 + i . 



2ivr(2 



(3.40) 



< — 



j f exp ((a + l/2)(i 2 - ln(i 2 ))) dt 
1 /"° 

y exp ((a + 1/2) (u 2 - 1 - ln(u 2 + 1))) dv = (3.41) 
exp(-( S + 1/2) ln(2)) ^ ^ ^ + m{ _ {i _ ^ + ^ _ i)2))) ^ ^ 

/•O 

y exp((s + l/2)(-r + 0(r 2 )))dr = (3.43) 

expf-r + 0(- — r — -)]dr= (3.44) 





2vr(2s + 1) 


exp(- 


-(s + l/2)ln(2)) 




2vr(2s + 1) 


exp(- 


-( S + l/2)ln(2)) 




?r(2s + l) 2 


exp(- 


-(s + l/2)ln(2)) 


tt(2s + l) 2 



.+1/2 (a + V2) 

2 



/ exp(-r)(l + 0( r ))dr= (3.45) 

A+l/2 (s + 1/2) 

= 0(exp(-(s + l/2)ln(2))). (3.46) 

Ainsi, sur les portions 7$ U 7_j du contour 7 , la contribution a l'asymptotique est expo- 
nent iellement petite : 

5SCSTI) eXP (<S + 1/2)( ' 2 " ln( ' 2))) dt = 0{ ^ h (3 ' 47) 

De meme, sur les portions (71Y71) U (7-A7-1) du contour 7 , la contribution a l'asymp- 
totique est exponentiellement petite par rapport a celle des portions limitees 71 U 7_i : 



> r-TXn / exp(( S + l/2)(t 2 -ln(t 2 )))dt 



(3.48) 



_ exp(s + l/2) ^ exp(-(2s + l) 2 / 6 ) , 

" ( 2s + l) 2 /3 ° ( (2s + W e ] ' ( } 

Ces deux ordres de grandeur etant exponentiellement petits par rapport a la contribution 
de l'integrale sur les portions 71 U 7_i , l'integrale sur le contour 7 vaut 



La methode du point col, pour la demonstration de cette proposition, permet aussi d'avoir 
un developpement asymptotique complet du coefficient [z 2s+1 ] T^(z) . Notons que dans 
les deux propositions, le terme principal est identique mais l'echelle asymptotique est 
differente : par l'expansion de Stirling, l'echelle est plus precise car c'est en 1/s, et par 
la methode du point col, l'echelle est, a priori, en 1/s 1 / 6 . 

La proposition 13.1.21 est un cas particulier du theoreme suivant : 
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Theoreme 3.1.3. Le nombre des hyperarbres (b-uniformes) enracines ayant (s(b — 1) + 1) 
sommets est determine a partir du coefficient de du developpement de la SGE 

T{z) = zexp(^^\ . (3.51) 
Ce coefficient est, pour (s — > oo) , tel que 



r s(b-l)+l 



n 



T(z) = i J -=\i + 0(4tz)} , (3.52) 



avec n = n(s) = s(b — 1) + 1 . 

Preuve. La definition implicite (I3.5ip donne, par la formule d'inversion de Lagrange avec 
la formule integrale de Cauchy, une expression integrale du coefficient |^z s ( b + 1 )+ 1 j de la 
SGE T(z) : 

avec n = n(s) = s(b — 1) + 1 et avec un contour qui encercle l'origine du plan complexe. 
Cette formulation integrale se simplifie en adoptant la notation r{t) que nous rappelons 
ci-apres 

r(t) = J^y / (3-54) 
En effet, nous obtenons l'expression suivante du coefficient 

[z n ] T(z) = (3.55) 

FT^qj / ex P (fzi) ex P W T W ~ ln ( r ^)))) T (3 ' 56) 



2iirn[(b 

1 



/exp exp (s(r(t) - ln(r(t)))) . (3.57) 



2wrn[(& - 2)!] s (6 - 1) J r ^ - V T W 
Comme la derivee logarithmique de la fonction r est 

£-<»-l)¥. (3.58) 

r definit un changement de variable, et pour garder un contour simple, il faut multiplier 
par (6—1) : 

^ T{z > = 2im [ (l !-2).]. / «■> (f^t) M'" T ' (3 ' 59) 

soit, en gardant la notation £ a la place de r , 
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Le parametre s dans l'exposant, nous indique que le point col a considerer est t = 1, 
soit la solution de l'equation 



±(t-ln(t)) = l-i=0. 
at t 



(3.61) 



Pour le choix du contour d'integration passant par ce point col t = 1 , nous prenons un 
demi cercle "gauche" 7 : 

71 :t = l + iv, v]e [-2,2] 



7 : 



70 : t = 1 + 2exp(ia) , a |e [tt/2, 3tt/2] . 



(3.62) 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 
-1 



1, v, 1-ln(1**2+v**2)/2 

1 +2*cos(u), 2*sin(u), 1 +2*cos(u)-ln((1 +2*cos(u))**2+(2*sin(u))**2)/2 




2 5 ^=^S^ 



2— +- 



FiG. 3.2 - La valeur de 3?(t — ln(t)) sur le contour (demi cercle gauche) 7 . 

Dans la figure I3.2I . nous voyons que la contribution , pour le coefficient recherche de 
la SGE, dans le module de Pintegrande sera tres accentuee au point col pour (s — > 00) . 
Notons k s tel que : 

ji:t = l + iv, v]e[-K s ,K s ]. (3.63) 
La contribution, au coefficient, de l'integrale sur cette portion de chemin de 7 est 



1 



exp 



2i7rn[(6-2)!] 
exp(s + l/(6- 1)) 
2inn[(b - 2)l] s ^ 



^-jj exp (s(t - \n(t))) y 



(3.64) 
(3.65) 



exp (s((t - l) 2 /2 + 0((t - l) 4 ))) dt , (3.66) 



soit, avec le changement de variable 



t = l + 



dt = ^dr. 



(3.67) 
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exp(s + 1/(6-1)) 



im[(b- 2)!] s 




exp(s + 1/(6 


-1)) 


irn[(b- 2)!] s 


y/2s 


exp(s + 1/(6 


-1)) 


7rn[(6- 2)!] s 


v 7 ^ 



exp 



-r 2 + 0(-; 



d r 



exp(-r 2 )(l + 0(— ))dr 
s 



exp(— r )dr + 



0(— )dr 



(3.68) 
(3.69) 
(3.70) 



Ainsi, sur la portion ji , approcher l'exposant donne Pordre de grandeur de Pintegrale : 
exp(s + 1/(6 - 1)) 



h 



im[(b- 2)\} s V2s 



exp(— r )dr + 0(- 



(3.71) 



En prenant k s = s 1//6 et en completant le domaine d'integration, pour obtenir la droite 
reelle dans Pintegrale et aboutir ainsi a Pintegrale de Gauss, il s'ensuit la contribution de 
Pintegrale restreinte a la portion 71 : 



«p(. + !/(>-!)) f 1+0( _L U . 
n[(6- 2)}} S V2^ \ V/ 6 ' r 



(3.72) 



Pour avoir le resultat, il faut encore negliger la contribution de Pintegrale sur le restant 
du contour defini en (|3.62p , a savoir 7i\tl puis 70 . 



1. Sur la portion 7i\tl du contour, la contribution au coefficient est 



1 



2»7ro[(&-2)!]- 
n exp( S + 1/(6-1)) 



cxp 

2s 



vrn[(6-2)!] s V2s 
= exp(s + 1/(6 - 1)) exp(-K? s ) . 
n[(6-2)!] s v / 2i k s 

2. Sur la portion 70 , considerons le chemin 



-M exp (s(t - ln(i))) y 



exp(-r 2 )(l + 0(— ))dr 
s 



(3.73) 
(3.74) 

(3.75) 
(3.76) 



7o : t = 1 + 2exp(ia) , a |G [n/2, 7r] , 



(3.77) 
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la contribution sur ce chemin au coefficient est 
h = 



2iirn[(b-2)\ 



J exp (j^Zj) ex P ( s (* - ln (*))) ~f 



<« ^-2W I »P«'-M*)))dt) 



o{ 



exp (l/(6 - 1)) 
i[(b-l 
e l/(6-l) 
vrn[(6-2)!] s 

e s(l+2cos(a)-ln((l+cos(a)) 2 +4sin(Q) 2 )/2)^ Q ,\ 



it/2 

e s-sln(VS)+l/(6-l) 



O f ^ ^ e-f^Dda) 



vrn[(6-2)! 



tt/2 



e 5-3ln(%/5)+l/(6-l) 

= ° ( vrn[(6-2)!] s ) = 

e s+l/(6-l) 

= °^im[(b-2)\\ a 5 a /^' 
Nous trouvons le meme ordre de grandeur sur la portion 7o\7o : 

Io = 

t \ dt 



/ ex P ( T~~- T ) exp - ln(i))) 



2inn[(b - 2)!]- J %Wq ~* \b-lj ~" ^ "™ t 

s -s H V5)+l/(b-l) rZn/2 

0{ TTr rrs / e 5 2, da 

1 7ro[(&-2)!]' h ' 

e s+l/(b-l) 

°^Trn[(b-2)l] s 5 s /^' 



3.78) 
3.79) 

3.80) 

3.81) 
3.82) 

3.83) 

3.84) 
3.85) 

3.86) 
3.87) 

3.88) 
3.89) 



Ainsi, la contribution au coefficient de la portion avec celle de la portion 70 sont 

exponentiellement petites par rapport a celle de la seule portion 71 du contour 7 . Nous 
deduisons, l'ordre de grandeur I\ du coefficient de la SGE : 



z s{b-i)+l 



exp I 5^ I ( 1 





Par definition d'une SGE A(z) enumerant des structures A etiquetees, le nombre de ces 
structures de taille n est n\ fois le coefficient [z n ] A{z) . Ainsi, connaissant l'equivalent 
asymptotique du factoriel par la formule de Stirling, n! ~ n n e~ n \/2irn , et connaissant 
un equivalent asymptotique du coefficient [z n ] T(z) par le theoreme 13.1.31 precedent, il se 
deduit le nombre asymptotique des hyperarbres enracines done de celui des hyperarbres 
non enracines (consequence directe de la proposition) . 
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Theoreme 3.1.4. Le nombre d'hyperarbres enracines ayant s hyperaretes et n = n(s) 
s(b — 1) + 1 sommets, pour (s — > oo), est : 



n\ [z n ] T(z) ~ 



Vb^Ts s ^ 



exp ((6 -2)( S + 1/(6-1))) 



1 



,6-1 



(6-2)! 



Preuve. Pour prouver ce theoreme il suffit de remarquer que 



ill 



n 



' exp(— n)V2nn , 



et que 



done 



[z n ]T(z) 



exp 



n 
6-1 



n! [z n ] T(z) 



n[(6-2)!]V 

/n/s 



2svr 



[(6-2)!] s exp (n 



n 
6-1 



(3.91) 
(3.92) 
(3.93) 

(3.94) 

❖ 



Et l'expression (|3,9ip s'ensuit. 

Le nombre des hyperarbres non enracines ayant n sommets differe d'un facteur n en moins 
que celui des hyperarbres enracines ayant le meme nombre de sommets, il devient clair 
que 

Theoreme 3.1.5. Le nombre d'hyperarbres ayant s hyperaretes etn = n(s) = s(b — 1) + 1 
sommets, asymptotiquement quand (s — > oo) ; est : 



n\ [z n ] o T(z) ~ 



,s(6-l)-l 



V6^Texp ((6 - 2)(s + 1/(6 - 1))) 



1) 



6-1 



(6-2)! 



(3.95) 



Ainsi, le nombre asymptotique des hyperarbres non enracines est determine indirecte- 
ment par la methode du point col car nous avons applique cette methode pour determiner 
l'asymptotique du coefficient [z n ] T{z) de la SGE relative aux hyperarbres enracines. Le 
point de depart de l'estimation asymptotique etant l'application de la formule d'inversion 
de Lagrange a la SGE pouvant etre exprimee en fonction de T(z) , pour avoir une expres- 
sion integrale du coefficient a estimer. Ayant a disposition la SGE des hypercycles lisses, 
l'estimation asymptotique du nombre des hypercycles est aussi accessible et elle est faite 
dans la section suivante. 



3.2 Enumeration asymptotique des hypercycles 

Dans cette section, pour enumerer asymptotiquement les hypercycles, nous estimons 
le n-ieme coefficient de la SGE Hq o T{z) que nous avons determinee dans le chapitre 
precedent (voir (|2.15p ) et qui est rappelee ci-apres 



H (t) 




t b-i 



f6-l 



2)! 



2)! 



(3.96) 
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Disposant de la formule d'inversion de Lagrange et de la formule integrate de Cauchy, par 
la methode du point col, nous determinons un equivalent asymptotique du coefficient de 
la SGE H o T{z) . 

Nous resumons la methode du point col, illustree par des preuves de la section 
precedente, comme suit : 

1. Fixer un contour passant par le point col (susceptible de convenir a la methode). 

2. Etablir le terme principal de l'asymptotique dans un voisinage, a preciser, du point 
col du contour. 

3. Negliger la contribution de l'integrale sur le chemin restant du contour, c'est a dire 
excluant la portion du voisinage defini plus haut. 

4. Conclure. 

Par une variante de cette methode, en autorisant que le contour ne passe pas exacte- 
ment sur le point col mais dans un voisinage proche, nous trouvons : 

Proposition 3.2.1. Le nombre des hypercycles ayant s(b — 1) sommets est determine a 
partir du coefficient de z^^ 1 ^ du developpement de la SGE Hq o T{z) avec 



Ce coefficient est, pour (s — ► oo) , tel que 

^ H ° ° T{Z) = ^(b-l^ {1 + 0{l/s)} , (3.98) 
avec n = n(s) = s(b — 1) . 

Preuve. La SGE Hq (|3.97|) des hypercycles lisses permet, grace a la formule d'inversion 
de Lagrange, d'exprimer le n-ieme (n = n(s) = s(b — 1)) coefficient [z n ] Hq o T(z) sous 
forme d'une integrate : 

avec un contour integral qui encercle, dans le sens direct, l'origine du plan complexe et 
r(t) = t^/ib - 2)1 . En effet, la derivee de la SGE H est 

Tt H ^ = wr^wr (3 - 100) 

Nous trouvons alors 

[z n ] H o T(z) = (3.101) 

1 / <*) . exp (p± - »MTM-m\ _ (3 . 102) 



MimJ l-r(t) ^V^- 1 b-1 J r(t) 

6—1 ft , , . . di 

exp(s(t-ln(t)))— = (3.103) 



4iirn[(b - 2)\} s J 1-t ^ y v wy/ t 
b-1 f 1 



4iirn[(b- 2)!] s J 1-t 



j> -^— t exp(s(t - ln(t)))d t . (3.104) 
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Sous cette forme integrate, nous observons que l'integrande presente un point col qui est 
determine par la racine de la derivee de l'exposant et un point singulier determine par la 
racine de (1 — t) . Soit t = 1 est a la fois un point col et un point singulier. Pour appliquer 
la methode du point col, nous adoptons la formulation integrate suivante du coefficient : 

l 2 "' "« ° T < 2 » " 4,4^-2)!]' / ^ (*<' " ' n(t> " ^f R) ) d ' ' (3 ' 105) 
Notons par h s le facteur du parametre s dans l'exposant de l'integrale precedente : 

h s (t) =t-\n(t)- kU-*) m ( 3-106 ) 

s 

Pour identifier d'eventuel point col, il faut trouver la racine de h s (t) = . La derivee h s ' 
etant : 

l "'^ = 1 -i + 7(r-r)' < 3 - 107) 

nous identifions le point col, de module strictement inferieur a 1 pour (s — > oo) , 



t< = 1 _vi?+i^i (3108) 

Nous faisons alors le choix d'un contour 7, similaire a celui de la preuve du theoreme 
13.1.31 . mais contournant le point col t s : 

r 71 : t = t s + exp(za)/2 , |€ [tt/2, 3tt/2] 
7: < 72 :t = t 8 + iv, v|€ [-3,-1/2] U [1/2, 3] (3.109) 

[ 70 : t = t s + 3exp(ia) , a |G [vr/2, 37r/2] . 

Nous voyons, sur l'exemple de la figure 13.31 . que le module de l'integrande capturera 
l'essentiel de Pasymptotique pour le parametre (s — ► 00) . Nous qualifions la methode 
comme methode du point col dans la mesure ou l'integration sur le contour 71 est egale 
a l'integration sur le segment t f € [t s — i/2, t s + i/2] . 

Le contour 7 d'integration (|3.109p etant fixe, il nous faut determiner le terme principal 
de Pasymptotique du coefficient porte par la portion 71 du contour autour du point col 
t s . Notons Ii la valeur de l'integrale sur la portion 71 , c'est a dire 

h = A b ~ l , n ( — e -^dt, (3.110) 

4ivrn[(6-2)!] s J^l-t t s ' v ' 

soit, avec le changement de variable 

dt^^ (3.111) 
71 ^ 71 : t = st s - s + se ia /2 , a |€ [tt/2, 3tt/2] , 
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j), 3*sin(u), .75+3*cos(u)-ln((.75+3*cos(u))**2+(3*sin(u))**2)/2-ln((.25-3*cos(u))**2+(3*sin(u))* 
5*sin(u), .75+.5*cos(u)-ln((.75+.5*cos(u))**2+(.5*sin(u))**2)/2-ln((.25-.5*cos(u))**2+(.5*sin(u))* 

.75, v, .75-ln((.75)**2+(v)**2)/2-ln((.25)**2+(v)* 
.75, -v, .75-ln((.75)**2+(v)**2)/2-ln((.25)**2+(v)* 


*2)/2/8 

*2)/2/8 

*2)/2/8 

*2)/2/8 



Fig. 3.3 - La valeur de 9?(i — ln(t) — ln(l — t)/s) sur le contour 7 pour s 

(b- l)exp(s) f 2 



4ivrn[(6- 2)!p 
(6 - l)exp(s) /" dt 



/ exp(0(t 2 /*))^ (3.112) 



4«7rn[(6- 2)!] 
(b — 1) exp(s 
4n[(6-2)!] s 

Ainsi, la contribution de l'integrale sur 71 est 



/ ^l(l + (t 2 /s)) (3.113) 
{1 + 0(1/ s)} . (3.114) 



II nous reste a negliger la contribution de l'integrale sur le contour restant, a savoir sur 
72 et sur 70 . Notons I2 l'integrale sur la portion 72 , alors 



6-1 



4inn[(b- 2)\} s 
(b — 1) exp(s) 
Amn[(b- 2)\} s 



(3.116) 

J exp(sh s (t))dt = (3.117) 
J 72 

/ exp (sh s (t) - s) dt . (3.118) 

J T2 



Rappelons l'equation parametrique de la portion 72 du contour ([3. l(J9j) : 

-f 2 :t = t s + iv, ute [-3,-1/2] U [1/2, 3]. (3.119) 
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Sur 72 , la partie reelle M(h s (t) — 1) , pour (s — > oo) , exprimee avec le parametre v de 
l'equation precedente, est telle que : d'un cote, si v € [—3, —1/2] , 

-1)<^(« + 3) -ln(VlO) (3.120) 
5 

et d'un autre cote, si v £ [1/2,3] , 

X(h a (t) - 1) < -^(« - 3) - ln(VlO) . (3.121) 
5 

La valeur de l'integrande li est done telle que 

h = (3.122) 

= q{ ( r7 1)e r p(s) /, r 1/2 ^ in(8)( ^ +3)/5 ^) + (3.123) 



v 47rn[(6-2)!] s 10 s / 2 7_ 3 

+Q( ( r ~ 1)e ^ (g) /9 T e'^^dv) = (3.124) 



v 47rra[(6- 2)!]nO s / 2 



1/2 



_ (6-l)exp(£) 

" Ul 47m[(6-2)!]-10»/2 J - ( ^ 5) 

Ainsi, la contribution au coefficient de l'integrale sur la portion 72 du contour 7 est 
exponentiellement plus petite que celle I\ sur la portion 71 . 

Nous trouvons aussi que la contribution de l'integrale sur la portion 70 , definie en 
(|3.1U9|) , est exponentiellement plus petite que I\ . Pour le montrer, notons Iq cette contri- 
bution : 

Jo = (3.126) 
/ exp (sh s (t)) dt = (3.127) 

■ / exp(sh s (t) - s)dt. (3.128) 

J Tn 



4t7rn[(6- 2)!]' 
(b — 1) exp(s) 



4iirn[(6-2)!] 



Rappelons l'equation parametrique du demi cercle 70 

7o : t = t 8 + 3exp(ia) , a |G [vr/2, 3vr/2] . (3.129) 

Avec cette notation, en prenant a comme variable, la partie reelle du facteur de s de 
l'exposant (h s (t) — 1) est telle que : 

- Dans un voisinage de t = t s + 3i , soit encore dans un voisinage de a = 7r/2 , 

9t(h.{t) - 1) = - ln(>/l0) - 2 l{a - 1) + 0((a - |) 2 ) . (3.130) 

- Et dans un voisinage de t = t s — 3i , soit encore dans un voisinage de a = 3tt/2 , 

R(^( t ) _ i) = _ in(VlO) + - (a - -) + 0((a - y ) 2 ) . (3.131) 
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Ainsi, l'integrale Io est telle que 

h = (3.132) 
_ (fr-l)exp(s) 

_CA 4™[(&-2)!]n(p/2 [6A66) 

I f e -f(-f)da+ r /2 e^ a -^da\) = (3.134) 

_ (6-l)exp(s) rt1 „,> 
" U( 4vrn[(6-2)!]n0^ j - 

Jo est done aussi exponentiellement petit par rapport a I\ . Nous concluons alors que I\ 
est l'ordre de grandeur du coefficient de la SGE et 

Wfl.°r M - (t i :' ) ^ (HO(:)}. (3-136) 



4n[(6-2)!] s [ V 



Notons que ce resultat est valide dans le cas des graphes en prenant simplement b = 2 et 
ainsi obtenir l'asymptotique du coefficient de la SGE Wo o T(z) des unicycles. 

Comme nous avons vu pour les hyperarbres enracines, le nombre asymptotique des 
hypercycles se deduit aussi de l'asymptotique du coefficient de sa SGE Hq o T(z) , donne 
par la proposition 13. 2. H . grace a la formule de Stirling : 

Theoreme 3.2.2. Le nombre d'hypercycles ayant s hyperaretes et n = n(s) = s(b — 1) 
sommets, pour (s — > oo), est : 



n\ [z n ] H Q o T{z) 



^(6-1) s^- 1 )- 1 / 2 
4 exp(s(6 - 2)) 



(6-1 



(6-2)! 



(3.137) 



Ainsi, nous avons une fois de plus la confirmation que l'analyse complexe, la methode 
du point col pour inspirer la recherche d'un contour d'integration, convient pour faire de 
remuneration asymptotique : par cette methode a ete etabli, theoreme 13.1.31 le nombre 
asymptotique des hyperarbres et , theoreme 13.2.21 . le nombre asymptotique des hyper- 
cycles quand la taille des structures est grande. Ces enumerations asymptotiques ont pu 
etre faites parce que nous disposons des SGEs des structures lisses correspondantes. La 
methode du point col, a priori, permettra aussi de trouver un contour aboutissant a un 
resultat d'enumeration asymptotique des composantes complexes d'exces £ > 1 donne, 
sachant que nous disposons d'un moyen "pratique" (un programme) pour determiner 1' ex- 
pression de la SGE Hi des structures lisses correspondantes. Dans la section suivante, nous 
procedons a remuneration asymptotique des composantes complexes selon leur exces. 



3.3 Enumeration asymptotique des composantes complexes 

Les preuves vues jusqu'ici pour avoir l'asymptotique des hyperarbres et des hyper- 
cycles necessitent la connaissance des SGEs des composantes lisses correspondant. Un 
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schema de preuve similaire permet aussi d'obtenir des resultats d'enumeration asymp- 
totique des composantes complexes, une des difficultes ici etant d'enoncer un resultat 
generique pour tout exces. 

Dans cette section, nous commengons par generaliser un resultat d'encadrement des 
coefficients des SGEs Wt o T{z) (cas des graphes) obtenu par Wright E.M dans [22| 
aux hypergraphes puis nous enongons un lemme pour l'asymptotique des coefficients de 
la SGE des sequences de m chaines, et enfin nous enongons le resultat d'enumeration 
asymptotique. 

3.3.1 Encadrement des coefficients de la SGE des composantes com- 
plexes 

Dans cette sous-section, nous tirons profit de la possibility, par le theoreme 12,3,371 
vu dans le chapitre 2, de determiner les premiers coefficients de la SGE, sous la forme 
(I2.120|) . des composantes complexes d'exces donne. Bien que dans ce chapitre, un langage 
analytique domine a des fins d'enumeration asymptotique, dans cette sous-section en 
particulier, nous retrouvons un raisonnement plus combinatoire (lecture bijective de SGE) 
comme dans le chapitre precedent. Nous soulignons en particulier, l'interet de la lecture 
combinatoire (voir la preuve du theoreme 12.3.21 ) de la formule d'inversion de Lagrange 
comme illustration cle pour avoir l'encadrement des coefficients. 

Pour les notations de l'enonce du theoreme qui suit, reprenons ici les notations utilisees 
pour les declinaisons des formes des SGEs Hi : 

H t (t) = (3.138) 

{i-e{t))n Jl fi-e(t) 



P =o v v ; 
= \i E Cj^b)0(ty= (3-140) 

(3.141) 

avec n = [(£ + l)/(b - 1) + lj et 9(t) = 1 - t 6_1 /(6 - 2)! . 

Dans le theoreme suivant, nous enongons l'encadrement des coefficients de la SGE 
Hi o T(z) des composantes complexes d'exces I > 1 . 

Theoreme 3.3.1. Le nombre n\ [z n ] HioT(z) de composantes complexes d'exces £ ayant 
n sommets admet l'encadrement suivant : 
La majoration 

n! [*»] Hi o T(z) < (n - 1)! \t^} { ^'^ Ht) n ) (3-142) 



et la minoration 



(« " !)> tn+l \ ( J^ETi ~ o^si ) ^ > * nl I*"] *l T ( z ) , (3-143) 
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ou 



et 



Or 



Bt = 3£{b - l)A £jSe , 
i 2 + i)A m + (b - 1)(3£ - l)(r e A ij3e - A 



e,3i-i) > 



= exp 



t 



(3.144) 
(3.145) 

(3.146) 



.En /ait, i?£ — C/> designent respectivement les coefficients de x~ 3i ~ l et de x~ 3i dans 

-(b-l)(l-x)f/(x)-£Mx), (3.147) 



3f 



1 — X 



(3.148) 



f e (x) = (l-xPY, A e 

p=0 

Preuve. Le nombre des composantes d'exces I > 1 ayant n sommets est determine grace 
a la formule d'inversion de Lagrange : 



n\ [z n ] H t o T(z) 



(n-1)!^ 1 ] 



$(i) n 



■(b-i)(i-e(t))ft'o6{t)-£f t oe(t)) 



(3.149) 
(3.150) 

(3.151) 
(3.152) 

(3.153) 

alors i?^ est un polynome de Laurent de degre minimum (—31 — 1) , de degre maximum 
borne par = \_(£+l)/(b— l) + lj , et a coefficients rationnels. Le nombre de composantes 
ci-dessus s'ecrit : 



= (n-1)! 
Notons pour I > 1 , 



($(t) n (-(6 - i)(i - 0(t))// o - eft o 



-(b-l)(l-x)//(x)-^(x) 



n! [z"]ff,oT(z) = (n-l)! 



(i? € o0(t)$(tn 



(3.154) 



II s'ensuit que les coefficients de la SGE R( o 9(f) sont tous positifs ou nuls, ce que nous 
notons pour une serie en la variable t comme suit : 



Rt°9(t) htO. 



(3.155) 



En effet, la multiplication par $(t) n qui apparait dans le second membre de (|3.154|) sou- 
ligne la decomposition des structures a enumerer, comme dans la preuve de la formule 
d'inversion de Lagrange (theoreme 12.3.21 ) : sur les structures lisses sont greffes des hy- 
per arbres enracines. Ainsi, Rn o Q(t) est une SGE enumerant des structures lisses IZi . 
Connaissant les caracteristiques du polynome de Laurent Rg , nous sommes assures de 
l'existence d'un nombre a^+i tel que 



0(f) 



(3.156) 
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En effet, il existe un entier N tel que les SGEs 

{9(tr 3e -\e(tr se , 9{t)-\ (1 - 9(t)), . . . , (1 - 9(t)) N } (3.157) 

forment une echelle de decomposition de Reo9(t) ^ . Combinatoirement, cette equation 
()3,156p dit que a%i + i9(t)~^~ l compte des structures lisses combinatoires beaucoup plus 
nombreuses a un nombre de sommets donne que TZi . Prenons a 3 i + \ le plus petit nombre 
verifiant cette propriete (j3.156H , formellement 

a 3e+1 = mm{a G 1R , ^l+i ~ R ? d ® ^ °> ■ ( 3 - 158 ) 

L 'equation (j3.156H souligne alors l'existence de structures combinatoires lisses (a une 
ponderation multiplicative pres) 1Z\ enumerees par la SGE R\ o 9(t) definie par son 
premier membre : 

R}(x) = ^ - Rt(x) . (3.159) 

Comme a^ + \ est choisi comme etant le plus petit, nous sommes assures de trouver un 
nombre tel que 

«'°m-(^-^f)t<0, (3.160) 

soit 

0L--Rjo9(t)htO. (3.161) 
Prenons ici aussi dans (I3.160P a 3 e comme etant le plus petit : 

a 3l = mm{a G 1R , - R\ o 9{t) y t 0} . (3.162) 

Le procede se reitere et avec un plus petit nombre a3^„i donne 
soit 

° 3/-1 R\ o 9(t) ht , (3.164) 



9(t) 3i ~ 1 



avec 



Rj( X ) = R e ( X ) - (ggl - gf) . (3.165) 

Le procede se reitere et s'arrete apres avoir determine a\ ou bien quand on obtient des 
structures TZ* 3 enumerees par la SGE Rf* o9(t) , Rl° etant un polynome. Ainsi, a un signe 
pres, Rj°(x) vaut 

*M-(|gi-|| + ... + (-ir- sf+1 |f). (3.166) 

En y faisant une identification des aj avec les coefficients du polynome de Laurent Ri , 
nous obtenons en particulier les deux inegalites du theoreme definissant l'encadrement 
comme consequence immediate de la propriete de positivite imposee a la definition des 
aj (qui par cette identification sont rationnels) : 
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la maj oration est, par (|3. 160|) , 

n! [z n ] H e o T(z) 
= (n- 1)! 

< (n- 1)! 

et la minoration est, par (|3.163p , 

n! [z rt ] H e o T(z) = 
= {n- 1)! [i n+l 

> (n-1)! r+n " ; 



f n+£ 



(Re o 0(i)*(t) 



< 



0(t) 



3£+l 



$(i) n ) 



o 0(t)$(*) r 



> 



(( 



Q31+1 
fl(t) 3 ^ 1 



Q3€ 



mt) n ) ■ 



(3.167) 
(3.168) 

(3.169) 

(3.170) 
(3.171) 

(3.172) 



Dans cette preuve, le raisonnement se base sur un procede d'inclusion exclusion pour 
compter des structures lisses qui contiennent un nombre fini de chaines : toutes les struc- 
tures sont comptees en exces comme si chacune avait 3£ chaines ; des extra structures 
seront alors indues dans le compte ainsi fait et ces extra structures peuvent a leur tour 
etre comptees en exces comme si chacune avait (31 — 1) chaines ; etc. 

Nous retiendrons la version de la majoration suivante : 

Theoreme 3.3.2. Le nombre n! [z n ] HioT(z) des composantes complexes d'exces i ayant 
n sommets admet la majoration suivante : 



n\ [z 11 ] H t o T(z) < 3lA m (b - l)(n - 1)! 



r(t) 



(l"T(t)) 



31+ 



(3.173) 



avec r(t) = t b ~ l / (b - 2)\ . 



Preuve. La preuve est similaire a celle du theoreme d'encadrement precedent en faisant la 
remarque que les structures comptees ont toutes au moins une hyperarete done la preuve 
est valide en faisant le remplacement Rg(x) <— Ri(x)/(l — x) dans l'equation (|3.153p . 



3.3.2 La contribution asymptotique de (m — 1) chaines 

Dans la justification du theoreme l3.3.2l . l'utilite de la classification des structures selon 
le nombre de chaines est mise en evidence pour proceder par inclusion exclusion et aboutir 
ainsi a la conclusion de l'importance, relatee par le theoreme precedent, de la contribution 
des structures a 3£ chaines dans le nombre des composantes d'exces t . La contribution 
dans l'asymptotique du coefficient de la SGE HgoT{z) provient essentiellement du terme 
avec le facteur 9(t)~ 3£ dans (|2.153|) . Cette remarque est intuitive par la forme meme de 
la SGE He et est appuyee bijectivement par le theoreme precedent. Cette relation, entre 
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les structures maximisant le nombre de chaines et l'asymptotique du coefficient de la 
SGE, a ete mise a profit par exemple dans [18] et dans [19] avec les notions de structures 
qualifiers de "clean" ou "unclean" . Pour tirer profit de la decomposabilite des structures 
en chaines, nous disposons du lemme suivant : 



Lemme 3.3.3. Pour (s — ► oo) , nous avons l'asymptotique du coefficient suivant 



un+t 



r(t) 



(1-T(t))» 



~2-b^(s/m) m l 2 



l + 0(- 



1 



,1/2 • 



(3.174) 



2^[(b-2)\Y 

avec n = n{s) = s(b - 1) - £ , r(t) = t b - 1 /(b-2)\, $(t) = exp(r(t)/(6 - 1)) et m > 2 . 
Preuve. La formule integrate de Cauchy donne 



r(t) 



1 



(l-r(t))" 
r(t) 



■ exp n 



fit) , tt1 „Adi 



2in J (I- r(t)) 
Comme n = n(s) = s(b — 1) — I , l'equation precedente devient 



Ltl+t 



r(t) 



(l-r(t))» 
1 



r(t) exp(sr(t) - sln(r(t))) dt 
2i7r[(6 — 2)!] s J (1 - r{t)) m exp^r(t) /(6 - 1)) T 



1 



r(t) exp (sr(t) - s ln(r(i))) d r(t) 
2ivr[(6 - 2)!] s (6 - 1) J (1 - r(t)) m exp(£r(t)/(6 - 1)) ^(t)~ 

exp (st — s ln(t)) 
2i7r[(6- 2)!] s J (1 - t) m " exp(£t/(6 - 1)) 



dt . 



(3.175) 
(3.176) 

(3.177) 
(3.178) 
(3.179) 
(3.180) 



ou le contour d'integration encercle, une fois dans le sens direct, l'origine du plan com- 
plexe. Avec l'integrale ainsi representee, nous remarquons de nouveau que 1 est a la fois 
point col et point singulier. Pour determiner un contour d'integration menant au terme 
principal de l'asymptotique, nous adoptons l'ecriture integrale suivante : 



r(t) 



(l-r(t)y 
1 

2z7r[(b — 2)!] 



-Hty 



exp (s(t - ln(t) - mln(l - t)/s)) 
exp(£t/(b- 1)) 



dt. 



(3.181) 
(3.182) 



Notons alors h , le facteur du parametre s dans l'exposant de cette expression integrale : 



hit) = t - ln(t) ln(l - t) . 

s 

Le point col a considerer est racine de h'(t) = . Comme la derivee h! est 



h'(t) = 1 



1 



+ 



t (s/m)(l-t) 



(3.183) 



(3.184) 
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le point col de plus petit module (strictement inferieur a 1) est 



vTOT-i . (3 , 85) 

2(s/m) 

Nous prenons alors un contour 7 similaire a (|3.62j) , celui qui a servi pour remuneration 
asymptotique des hyperarbres enracines : 

f 71 :t = t + iv/s, v]£[-3s,3s] ^1861 

\ 7o :t = t + 3exp(ia), a |G [tt/2, 3vr/2] . l " J 

Soit une portion 71 du chemin 71 : 

71 :t = t Q + iv/s, v ]£ [-k s ,k s ] , (3.187) 

avec n s , un nombre positif qui sera precise plus tard. Notons I\ la valeur de l'integrale 
sur cette portion 71 : 

71 -2^(6-2)!]-;^ exp(£t/(6-l)) dt ' (3 " 188) 
soit, avec le changement de variable 

dt^ ^ (3.189) 

7i l— > 7i : * = sto — s + ^ j v |G [—3s, 3s] , 



exp(-&) 



avec 



Comme 



/ S exp(s/(u))du, (3.190) 

J— K. 



2stt[(6-2)!] 

/(«) = /i(to + w/s) - £iv/(bs 2 - s 2 ) . (3.191) 



»(/(«)) = (3-192) 
= ft(/i(to + iv/*)) = (3.193) 

= to - i ln(tg + (^/s) 2 ) - ™ ln((l - t ) 2 + (v/s) 2 ) = (3.194) 

m , , s . m „. 1 . , 
= 1 + - In - + - + O {— j + 3.195 

2s m 2s s d / z 

u 2 , 3m 1//2 5m lira 3 / 2 . 1 . . 

-?< 1 + W + V + Tii5^ + 0( ? ))+ < 3 ' 196 > 

+0(^j), (3.197) 
s 6 
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h = 



„ , m l*Q 

e + 2 fc- 1 



s/m) 



m/2 



2sir[{b-2)\ 

x 



K s u 2 /i i 3m 1 / 2 | 5m , 11m 3 / 2 

2 S V 2 



+ ^7 2 - ) (l+0(-^))dt;. 



Notons 



1 , 3m 1 / 2 5m llm 3 / 2 



16s 3 / 2 



alors par le changement de variable 



v = 



/C(s) 



C(s) ' 



l'integrale I\ est telle que 
h = 



e s+ 2 


— ( s / m )m/2 


2stt[(6 


-2)!]V^) 


e + 2 


^(s/m) m / 2 


2stt[(6 


-2)!]V^) 


e + 2 


^1 (s/m) m l 2 


2svr[(6 





_e~ r '(l + 0(_±_))dr 



e~ r dr + 0( 



C(s)s 3 ' 

- 3 C(s) 3 / 2 , 



C(s)s 3 



e r dr + 0(— Y — ) 



Comme pour (s — > oo) , C(s) ~ 1/s , nous obtenons 



2stt[(6-2)!] 



Ks/t/s 



e~ r dr + 0(- 



: 7/2 ' 



(3.198) 
(3.199) 

(3.200) 
(3.201) 



(3.202) 
(3.203) 

(3.204) 

(3.205) 



(3.206) 



Cette ordre asymptotique est valide, pour m > 1 fixe, en prenant k s = s et se reecrit 
alors : 



_ e s+ ^-^(s/m) m / 2 ^ j f 
1_ " 2stt[(6-2)!]- \J_ 



e~ r dr + 0(- 



1 



,1/2 - 



(3.207) 



Ainsi, Pequivalent asymptotique de l'integrale I\ sur le segment 71 de longueur deux fois 
(k s = s) est : 



h = 



e s +f-^( s / m )W2 



l + 0(- 



1 



,1/2 < 



(3.208) 



2^(6-2)!]* 

Pour obtenir le resultat, il faut encore negliger la contribution de l'integrale sur le contour 
rest ant. 
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1. Sur la portion 71Y71 , la contribution I[ est 

j' _ 1 / exp(sfr(t)) 

avec 

7i\7i : t = t + w/s , v |€ [-3s, — /c s ] U [k s ,3s] (3.210) 
ou bien (comme nous avons pris k s = s ) 

7i\7i :* = *o + *u, ut€ [-3,-1] U [1,3]. (3.211) 
L'integrale /[ , pour (s — > 00) , devient alors 

s+m/2-Q, , >. m /2 , ,-1 ,3^ 

'! = CT^F {/., + / I exp(s/( " ))d " ' (3 ' 212) 

avec la fonction / definie comme suit 

t, . . ^ m(l + ln(s/m)) ln(s) 

/W = fc(ta+w , ) -_- 1 - _JJ (3.213) 

Sur l'intervalle i> G [—3,-1], nous pouvons caracteriser la partie reelle $t.(f(v)) 
comme suit : 



= R(Afo + it,) - 1 - m(1 "^ /m)) + ^ < (3.215) 

2s 2s 



»(/(«)) = (3.214) 
R(/i(to + iv) - 1 - - 

<!2S) ( „ +1) _!^. (3 . 216 ) 

Et sur l'intervalle t> G [1,3] , nous obtenons la caracterisation suivante : 

»(/(«)) = (3.217) 

= + *,)-!- m(1 ^f /m)) + ^ < (3.218) 

<-^-D-^. S S (3-219) 
Des deux majorations precedentes, nous en deduisons que l'integrale I[ verifie 

I[ = 0{- F ptt^S (3.220) 

1 v 2 v /i 7 r[(6- 2)!] s 2 s / 2 v ' 

-1 /-3 



jy e .ln(2)(H-l)/4 dt; + jf e - a ln(2)(«-l)/4 du |) (g^l) 



et l'ordre de grandeur asymptotique de l'integrale J( est 

s+m/2( / \m/2 i 

^ VSW 2).]. ^' <3 ' 222) 

c'est a dire que l'integrale /J est d'un facteur exponentiellement negligeable par 
rapport a I\ . Nous montrons de meme pour le contour restant. 
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2. Sur la portion 70 , definie a Pequation (|3.186p , la contribution Iq est 



1 



exp(sh(t)) 
2»tt[(6-2)!]- 7 70 exp(W-l)) 

z s+m/2 -^{s/m) m / 2 



dt 



2^[(b -2)!] 



+ 

tt/2 ./tt 



3tt/2 



exp(s/(a))d a . 



avec la fonction / definie comme suit 



f{a) = h(t + 3e i 



3£e ia 

+ Z--1 

so — s s 



m(l+ln(s/m)) m (s) 



+ 



2s 2s 

Sur l'intervalle a € [7r/2, tt] , nous avons la partie reelle 3ft(/(a)) telle que 

»(/(«)) = 

= »(M*o + 3e'«)) - - 1 - !^l±Mi» + K-) 



s6 — s 



2s 



2s 



= - ln(VlO) - + C(-^))(a - |) + O(-L) + 0((a - ^) 2 ) . 

Et sur l'intervalle a £ [7r, 37r/2] , nous avons la partie reelle 5R(/(a)) telle que 
»(/(«)) = 

= - ln(>/lO) + + 0(i=))(a - + O(-L) + 0((a - ^) 2 ) 
Aussi, l'integrale Io est telle que 
io = 

gS+m /2-«!L (s/m)m/2 j 



2yivr[(6-2)!] s 10 s / 2 

e -27 S (^/2)/io da + 

tt/2 

s+ m /2-^ (s/m)m /2 j 



f 



3?r/2 



,27s(a-37r/2)/10 d 



a 



2 v /stt[(6- 2)!] s 10 s / 2 



) 



Jo est done aussi exponentiellement petit par rapport a Ji . Nous concluons que I\ est 
l'ordre de grandeur du coefficient : 



r(t) 



(l-r(t))' 



-d>(i) n 



1 m __^0_ 

2v^F[(6-2)!] 



to/2 



i + o(- 



1 



,1/2 • 



(3.236) 

❖ 



Ce lemme offre la possibilite de faire le saut entre remuneration exacte de nos structures 
decomposables en chaines et leur enumeration asymptotique. 



62 



CHAPITRE 3. ENUMERATION ASYMPTOTIQUE 



3.3.3 Enonce du theoreme d'enumeration asymptotique 

Proposition 3.3.4. Le coefficient [z n ] HioT{z) du SGE des composantes d'exces I admet 
V equivalent asymptotique suivant 

[z n ] H e o T(z) = (3.237) 
avec n = n(s) = s(b — 1) — I , ou s est le nombre d'hyperaretes. 

Preuve. Par le theoreme 13.3.21 . l'equivalent asymptotique du coefficient recherche est porte 
par 



I = 3lA m ^ 
n 



t n+t 



(l-r(O) 



avec r(t) = t b - l /(b - 2)! et = exp(r(t)/(6 - 1)). Et par le lemme l3~3~3l . avec 

m = 3£ + 1 , nous obtenons l'equivalent asymptotique de I suivant : 



j | 31+1 £ , a , 31+1 







Theoreme 3.3.5. Le nombre n\[z n \Hi o T{z) des composantes complexes d'exces I, 
ayant s hyperaretes et n = n(s) = s(b — 1) — I sommets, admet l'equivalent asymptotique 
suivant : 



n\ [z n ] H £ o T(z) = (3.241) 



/2^ e s6- 2s +^/(6-i) [(6-2)!]* [ V V« 



-A^3^ defini dans la forme (I2.22D de Za 5"Gi? , se determine par le theoreme 



2.3.31 
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Preuve. Nous avons par le theoreme 13.3.21 et par le lemme 13.3.31 avec m = S£ + 1 
n\ [z n ] H e o T(z) = 

s | se+i e_ 3t+i 

= 3( A tMb -!)(„- l / {l + °<£) 



es -zt+i 



,ra N „ e ^ 3^+1 J I - 1 



3M £ 3£ (6 - l)(-) n ; ^ 1 + 0(-=) 



^< es s^+i 



= uA t ,^ - i)(-r e ~^ +1 LI {1 + o(±: 

e V2sn[(6 - 2)!] s [ V s 

-^O-D^- L^!. { 1 + 0( -L; 

^ e (p-2)/(b-i)> ^/2i^[(&-2)!]s I V ^F 

Comme pour (s — > 00) 

n n = (s(b - 1) - £)< 6 - 1 )^ = [a(6 - l)]^ 6 - 1 )-^ exp(-l) , 
l'equivalent asymptotique recherche s'ecrit 
n\ [z n ] H e o T(z) = 

-MA (b n ML^^ ( ^ } ^ fury 1 

31+1 „ N „ 

1 



2sne sfe-2 S +£/(fe-l) [(6-2)!] s [ 



<> 



Dans ce chapitre, qui est l'enchainement de l'enumeration exacte, des enumerations 
asymptotiques ont ete obtenues avec un recours a Panalyse complexe. Nous soulignons en 
particulier Putilisation de la methode du point col pour toutes les enumerations asympto- 
tiques faites : des hyperarbres enracines aux hypercycles et aux composantes complexes 
selon leurs exces. Notons que le chemin que nous avons choisi ici devoile completement 
la "magie" de la formule de transfert donnant a partir des SGEs et de leur forme, le 
terme asymptotique principal voire l'expansion complete de ses coefficients : en parti- 
culier, les contours d'integration choisis peuvent servir pour expliciter des expansions 
completes. En perspective immediate du travail fait dans ce chapitre est la consideration 
de l'enumeration asymptotique des composantes complexes d'exces infini et suffisamment 
peu denses pour que les preuves donnees ici restent valides - un tel resultat parait dans 



64 



CHAPITRE 3. ENUMERATION ASYMPTOTIQUE 



[21, 27J. Ainsi. la plupart des travaux faits sur lcs graphcs pent trouver une generalisation 
par une demarche similaire a celle suivie jusqu'ici : passant par une recurrence des SGEs, 
identifiant la forme des SGEs, puis par la formule d'inversion de Lagrange et la formule 
integrale de Cauchy, etablir une expression integrale des coefficients et en deduire alors 
une portion du contour capturant la contribution principale de l'integrale. 

Dans le chapitre suivant, nous illustrons par les hypergraphes aleatoires qu'effecti- 
vement, les travaux et resultats existant sur les graphes peuvent etre generalises aux 
hypergraphes et ainsi, revisites et compris via les SGEs et l'analyse complexe pour obte- 
nir des caracteristiques asymptotiques des structures. 



Chapitre 4 

Hypergraphes aleatoires 



"On s'appuie de l'histoire. Mais notre histoire n'est pas notre code. Nous devons nous 
defier de prouver ce qui doit se faire par ce qui s'est fait. Car c'est precisement de ce qui 
s'est fait que nous nous plaignons" (1788). 



Dans ce chapitre, nous adoptons une vision dynamique des hypergraphes : dans la 
prochaine section, les hyperaretes d'une composante sont recursivement enlevees jusqu'a 
ce qu'il n'en reste plus et dans une autre section, les hyperaretes sont ajoutees une par 
une jusqu'a l'obtention d'une certaine propriete dans la structure. Notre outil principal 
pour mener notre etude restent les SGEs (voir [301 [2.6] ) qui permettent par exemple dans 
[T3] d'obtenir diverses caracterisations statistiques. 

4.1 Hypergraphes et hypercouplage glouton 

Cette section est une generalisation aux hypergraphes de resultats de [H] sur l'analyse 
en moyenne de la performance de l'algorithme glouton de couplage. 

Definition 4.1.1. Un hypercouplage est une collection d'hyperaretes deux a deux dis- 
jointes. 

Le probleme de trouver dans un graphe un couplage qui maximise le nombre de ses 
aretes est connu en informatique. Ce probleme admet une version avec les hypergraphes : 
maximum hypercouplage, consistant a trouver un hypercouplage, ayant un nombre maxi- 
mum d'hyperaretes, dans un hypergraphe. 

L'algorithme [3] est un algorithme glouton et il fournit un resultat qui, de maniere 
generale, n'est pas optimum, mais presente l'avantage de la rapidite et garantit tout de 
meme que 1' hypercouplage qu'il retourne soit au moins maximal a defaut d'etre maximum 
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Algorithme 3 : Hypercouplage glouton sur hypergraphe 
Entrees : Une composante. 
Sorties : Un hypercouplage maximal, 
debut 

Initialiser avec un hypercouplage vide, 
repeter 

Choisir aleatoirement une hyperarete a placer dans Phypercouplage. 

Supprimer cette hyperarete et celles qui lui etaient adjacentes. 
jusqu'a il n'y a plus d'hyperarete 
retourner Vhypercouplage construit 

fin 



pour l'hypergraphe. Notons aussi que c'est un algorithme non deterministe car le choix 
d'une hyperarete a ajouter dans Phypercouplage est aleatoire. 

Dans cette section, nous proposons de faire l'analyse de la performance de cet algo- 
rithme glouton en determinant la taille moyenne des hypercouplages qu'il retourne, soit 
le nombre moyen d'hyperaretes contenus dans ces hypercouplages. Une telle analyse a 
ete faite dans pj] pour le cas des graphes. Dans la sous-section qui suit, est presente le 
formalisme mathematique pour proceder a cette analyse de l'algorithme glouton. 

4.1.1 Definitions et notions 

Dans cette section, notre but est d'apprecier la variable aleatoire qu'est le nombre 
d'hyperaretes contenues dans Phypercouplage produit par l'algorithme [3] glouton, quand 
l'hypergraphe de depart est une composante d'exces I et ayant un grand nombre s d'hy- 
peraretes (done un grand nombre de sommets n = n(s) = s(b — 1) — I) . Nous fixons ici 
les notations pour faire l'analyse : 

(£) 

- Yn ■ la variable aleatoire correspondant a la taille de Phypercouplage produit par 
l'algorithme glouton quand Pentree de l'algorithme est choisie uniformement parmi 
les composantes d'exces i , ayant n sommets et s hyperaretes tels que n = n{s) = 
s(b-l)-£. 

- fn^ '■ la fonction generatrice de probabilite (FGP) associee a la variable aleatoire 

1 n j 

fW(z) = Y,P(YW=y)zy, (4.1) 
y>o 

ou 

- z denote la variable liee au nombre d'hyperaretes d'un hypercouplage. 

La derivee de la fonction generatrice de probabilite (|4.ip evaluee en z = 1 est Pesperance 

<£) , <£) 

En de la variable aleatoire Yn , soit du nombre d'hyperaretes dans un hypercouplage 
produit par l'algorithme glouton quand Pentree est choisie uniformement parmi les com- 
posantes d'exces £ ayant n sommets. Comme la fonction de probabilite vaut 1 quand 
elle est evaluee en z = 1, pour etudier la recurrence sur n arm d'obtenir une estimation 
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Fig. 4.1 - Une trace d'un deroulement de Falgorithme [3] sur une composante explicitee 
dans l'annexe. 



asymptotique de l'esperance Yn , nous introduisons la SGE G^ bivariee 

GW(x,,) = ^ c W/W(.): 



n>0 



??■! 



(4.2) 



avcc 



Sa '■ le nombre de composantes d'exces I ayant n sommets, 

x : denote la variable liee au nombre de sommets de l'hypergraphe d'exces I donne 
en entree. 
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Une propriete des FGPs fn nous amene a definir la fonction generatrice de la moyenne : 



EW{x) = ^-G®{x,z) 



1 n: 



Nous avons aussi la propriete de la SGE G(£) suivante : 

G®(x,l)=H e oT(x), (4.4) 

avec Hi o T(x) , la SGE des composantes d'exces £, et T(x) , la SGE des hyperarbres 
enracines, dont nous rappelons la definition implicite comme suit : 

r(x)=xexp (^i9) ■ (4 - 5) 

Remarque 4.1.2. Notons la traduction, soit la lecture par rapport aux deux variables, 
des operations sur les SGEs (|4.2j) bivariees dans 

C&(x,z) + gM(x,z) = £teW/„ W (*) + c^fi k \z))- (4.6) 

n 

G®(x,z) x G W(x,z) = ( f\c?c^\z)fi%(z)^ . (4.7) 

La lecture combinatoire pour les sommets etiquetes pour les deux operations d'addition 
et de multiplication suit la justification du dictionnaire de |26| . Une lecture combina- 
toire comme l'union et le produit, est alors immediate par les proprietes d'une fonction 
generatrice de probabilite si les variables considerees sont independantes. Aussi, pour uti- 
liser les SGEs bivariees G^' , il est important de garantir cette independance des variables 
aleatoires considerees dans la decomposition souhaitee. 

Dans la prochaine sous-section, nous donnons la decomposition qui nous servira pour 
1' analyse. 



4.1.2 Decomposition ou recurrence 

Dans le chapitre sur remuneration exacte, nous avons etabli une recurrence des SGEs 
Hi des composantes d'exces I en partant de la decomposition suggeree par le marquage 
d'une hyperarete : c'est la decomposition utilisee par Wright dans [22J pour les graphes. 
Cette decomposition ne se traduit pas facilement (directement) avec les SGEs bivariees 
car une fois la decomposition faite, nous perdons la lecture combinatoire dans les com- 
posantes dissociees : les variables aleatoires ne correspondent pas aux differentes tailles 
des hypercouplages produits avec ces composantes en entree de Palgorithme glouton. 

L'algorithme[3]d'hypercouplage glouton suggere une nouvelle decomposition des struc- 
tures : une composante d'exces I ayant une hyperarete marquee, se decompose, non pas 
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seulement en ignorant Phyperarete mais aussi toutes les hyperaretes qui lui etaient ad- 
jacentes, en des composantes d'exces plus petit. Pour etablir une bijection, il nous faut 
pouvoir recombiner les composantes dissociees, en la composante d'exces £ de depart. 
Pour proceder bijectivement a cette recombinaison des composantes dissociees, nous de- 
vons marquer un ou plusieurs sommets de ces composantes pour les lier a Phyperarete 
marquee en creant d'autres hyperaretes, chacune pouvant avoir bo = 1, . . . , b — 1 som- 
mets en commun avec cette hyperarete marquee, k = 1, . . . , b — bo sommets (parmi ceux 
marques) en commun avec une ou plusieurs composantes dissociees et b — bo — k racines 
d'hyperarbres. 

La decomposition est done descriptible par des ensembles de composantes ayant k = 
1, ... ,b — 1 sommets marques. Par exemple, dans le cas ou l'exces vaut I = — 1 , soit 
Panalyse de Palgorithme quand en entree un hyperarbre est choisi uniformement, une 
telle decomposition se traduit par la proposition qui suit : 



Proposition 4.1.3. Le marquage d'une hyperarete d'un hyperarbre se traduit comme 
suit .' 

JL r(l) + ^ TW ) _ £ exp („Zg_' 

Proposition 4.1.4. Pour tout £ > , la SGE Hi des composantes d'exces £ satisfait la 
relation : 

£H e o T(x) + x^-Ht o T(x)) = ^ exp ( b^—) x (4.9) 



6-1 V — K ' dx v 7 6! "V (6-1)! 



£ b-1 k 



Cyc e+1 \ exp I ^^ C yc^ fc /3A ; (Cyc,x)x fc ^^or(x)+ (4.10) 

j=0k=l 



b ~ l i k 



+ Y, Cyc^ 1 A(Cyc, x )x k ^H_ l o T(x) ) , (4.11) 

k=2 



avec 



p k (Cyc,x)= (4.12) 
1 fT(x) b - 2 \ k A 1 / u T(x) b -* V 

»X b (^w) ^ Q A pb ^ Cyc ) + (413) 

\b-2\ k - 1 /Tf„\b-3\ e_1 1 / Tf~\b-2\P 



p=0 

2 J (k — 2)! 

+ ... (4.16) 

Les decompositions donnees dans ces deux propositions sont valides en notant que 
T(x) = -^H-i o T(x) , puis en substituant les o T(x) par les SGEs G^(x, z) corres- 
pondantes (des derivees droites en x deviennent des derivees partielles) a un facteur z 
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pres car une hyperarete a ete choisie et mise dans l'hypercouplage : la decomposition est 
en accord avec l'independance des variables aleatoires dans les differentes composantes . 
Nous avons done : 

Proposition 4.1.5. Dans le cas des hyperarbres, 



h l 



(- G H)(v) + 4gH)(v)) =4^( & ( ^V-iT )) " • (417) 



Entre autres motivations, afin d'alleger les notations, nous adoptons les notations des 
operateurs de marquages suivantes : 

Definition 4.1.6. Dans des structures, nous traduisons le marquage de k sommets or- 
donnes par Poperateur sur des SGEs 

( 4 - i§ ) 

si x designe la variable liee aux sommets. 

En particulier, pour le marquage d'un sommet 

(419) 

Et nous adopterons cette notation de marquage selon la notation de la variable utilisee : 
tit = tm ou encore ti y = y-j^ . 

Definition 4.1.7. Dans le cas des structures d'exces I, nous traduisons le marquage 
d'une hyperarete par Poperateur sur des SGEs 

*« W = ^ , (4.20) 
si w designe la variable liee aux hyperaretes et x , celle liee aux sommets. 

4.1.3 La performance gloutonne sur les hyperarbres 
Theoreme 4.1.8. La SGE 

9-1 (y) 



E^ 1 \x)=T{2 



(l+y)l/(6-l) 



y (6-2)! 



(4.21) 



de la moyenne de la variable aleatoire Yn est determinee par la recurrence 



(6-l)(l + y)(f-2)/(&-i) 5 _ 1 '(y) = \, (4.22) 
par laquelle g~i(y) est uniquement determinee avec la condition <7-i(0) = . 
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Preuve. En differenciant par rapport a la variable z, l'equation de la proposition prece- 
dente, puis en fixant la valeur z = 1 , nous obtenons 



2_(_^ G <- 1 ) (I , Z)+ , I | G (- 1)(X ,, ) ) 



Z = l 



b\ 



exp b- 



(6-1)! 



+z b 



(^ x G^(x,z)) b - 
(6-2)! 

x b 

x — exp ( 6 



(6-1)! 



2=1 



soit 



tiJ-VE^ix) 
— exp o 



6! 

T(x) b 
6! 



Tjxf- 1 
'(6-1)! 



1+ 6 



T(x) 



1+ 6 



T(x 



,6-2 



(6-2)! 



(6-2)! 



soit en passant les occurrences de E(~V au premier membre 



b(-E^~ 1 \x) + (l 



T(x 



,26-2 



[(6-2)!p 



VxE^Hx) 



T(x) b 
(6 = 2)! 



ct 



T(x) 



+ 1 



T(x) 



6-1 \ 2 



(6-2)! 



^(^(x) _ 1 r(x) b ~ 1 
T(aO ~ 6 (6-2)! 



Determinons alors ^ (x) sous la forme 



E^ 1 \x)=T( 



x 



g-i(y) 



(l +y )l/(6-l) 



" (6-2)! 



Notons que pour une fonction / , si r(t) = t /(6 — 2)! , 



V x (T(x)/oroT(x)) = 

= (^T(x)) x (/ o r o T(x) + T(x) x r' o T(x) x /' o r o T(x)) 



T(x) 



y=roT(x) 



(4.23) 
(4.24) 

(4.25) 

(4.26) 



(4.27) 
(4.28) 

(4.29) 



(4.30) 



(4.31) 



(4.32) 



(4.33) 
(4.34) 

(4.35) 
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(I4,32[) permet le changement de variable 
( T{xf- 1 



g-i(y) 



(6-2)! 

T{x) ~* (l + y) 1 /^-!) 
•& X E(-V{x) , /l + (6-l)^\ 



(4.36) 



T{x) V 1-2/ y (l+y) 1 /^" 1 ) 

qui, si nous le portons dans Pequation (|4.3ip . donne 



(l + y)V(6 
D'ou nous deduisons 



g-i(y) 



(i + y )i/(6-i) y 6 



(6-l)(l + ^- 2 )/^ 1 ) 5 _ 1 '( y ) = i. 



(4.37) 
(4.38) 



La resolution de la recurrence, que definit Pequation differentielle (I4.22H , est immediate 
et nous donne : 

Theoreme 4.1.9. Dans le cas des hyperarbres, V expression de la SGE, de la moyenne 
de la variable aleatoire Yn ^ , est 

T(x) 



E^(a 



1 



V 



1 + 



T(x) b 



n\ V(6-l) 



(6-2)1 



(4.39) 



Notons que ce resultat generalise celui de [IT] qui est le cas des arbres (6 = 2), dans 
quel cas, 1' expression de la SGE est 



£(-!)( 



T(xf 



2(1 + T{x)) 



(4.40) 



Disposons de la SGE E^ 1 \x) , de la moyenne de la variable aleatoire Yn ^ , qui s'ex- 
prime en la serie T(x) des hyperarbres enracines, nous sommes en mesure de proceder a 
l'analyse du comportement asymptotique de cette moyenne via le comportement asymp- 
totique du coefficient de la serie. 

Theoreme 4.1.10. L 'equivalent asymptotique du coefficient [x n ] E^^(x) de la SGE liee 
d la moyenne de la variable aleatoire Yn ^ est : 

, e n/(b-i) (1 _ b/2 b/(b-i) ) 



nb[{b-2)\} s V2TTs 



l + 0{ 



,1/6 - 



(4.41) 



avec s le nombre d'hyperaretes et n = n(s) = s(b — 1) + 1 . 
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Preuve. Notons E*> 1 ) (t) la SGE lisse correspondant a (x) . Alors, E^ ^ (x) = & ^ o 
T(x) et 

/ 



E^\t) 



La derivee de la SGE E^it) est 



1 + 



(6-2)1 



1/(6-1) 



dt 



E^Xt) 



1 



+ 



r(t) 



b bil + rit)) 11 ^ b{l + r{t)) b/{b - 1) 
1 1 



b 6(l + r(t)) 



6/(6-1) ' 



avec r(t) = t b ~ l /(b - 2)! . Notons alors 



(4.42) 



(4.43) 
(4.44) 



E^\y) 
1 



1 



1 1 



b 2 b /Q>- 1 )bJ (b- l)2 1 /(6-i)2 2 



T^(y-i) + - 



(4.45) 
(4.46) 

(4.47) 



Soient s , l'entier correspondant au nombre d'hyperaretes et n = n(s) = s(b — 1) + 1 , le 
nombre de sommets, alors la formule d'inversion de Lagrange donne le coefficient 



[z n ]E^ l \x) = 

1 r p exp(ftr(t)/(6 - 1)) dt 



2%-nn 



1 



2t7rn[(6-2)!] 



M-V or(t) 



t n ~ l t 

exp(nr(t)/(6 - l))drft) 
(&-l)r(t)« r(i) ' 



(4.48) 
(4.49) 

(4.50) 



avec un contour d'integration, suffisamment petit, qui encercle une fois l'origine dans le 
sens direct. Ainsi, 



[x n ] E^ia 



1 



2«rn[(&- 2)!] 
1 

2iirn[(b- 2)!] 



|,(-i)^ exp(nt/(6 - 1)) df 



^(tJexpCt/Cft-l)) 



exp(si) dt 

i* r 



(4.51) 
(4.52) 

(4.53) 



L'expression de la derivee (|4.45p nous indique que l'ordre de grandeur asymptotique 
recherche est celui de l'integrale 



1 - l/(2 fo /(f'-i)) r 
~ 2nbi7r[(b - 2)[] s J 



ex P (V(6-l))^^di, 



(4.54) 
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qui correspond au facteur pres a (|3.60p . Aussi nous obtenons 

nb[(b- 2)\} S V2^ X V/ 6 7 y ' 

soit l'ordre de grandeur recherche. 



Connaissant le comportement asymptotique du coefficient de la SGE E^~ 1 \x) de la 
moyenne de la variable aleatoire Yn ^ et celui de la SGE H^\oT{x) , deduite du theoreme 
13.1.31 . par la definition (|4.3p , il decoule : 

Theoreme 4.1.11. La moyenne En ^ de la variable aleatoire Yn ^ correspondant a la 
tattle de I'hypercouplage produit par Valgorithme glouton quand I 'entree de I'algorithme 
est choisie uniformement parmi les hyperarbres ayant n sommets, est telle que 

^-(l-l/^ 1 ))^. (4.56) 

4.1.4 La performance gloutonne sur les £-composantes (£ > fixe) 

Introduisons, pour l'analyse dans le cas des structures d'exces £ > en particulier, les 
fonctions F^(x, z) , pour reduire la longueur des formules decrivant des decompositions. 

Definition 4.1.12. Notons la SGE bivariee 

F^(x,z) = (4.57) 

(l 6-1 
Yj Q y cj+k Pk ° (x, z)$ x k G^ (x,z)+ (4.58) 
3=0 k=l 

6-1 \ 

+ Y,Cyc k - 1 (3 k oG(- 1 \x,z)$ x k G(- 1 \x,z)\ , (4.59) 

k=2 J 

avec 



k (Z) = (4.60) 



« V (b-2)\J f^pl \T (6-2)! 
+ ... (4.64) 
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Cette definition des SGEs permet de formuler de maniere concise les equations 
des propositions 14,1.31 et 14.1.41 en evaluant l'equation suivante en z = 1 : 



x b fj x G ( -- 1 \x 1 z) b - 1 



b\ 



F®{x,z), 



(4.65) 



(6-1)! y 

qui est l'ecriture des decompositions relatees dans ces propositions avec les SGEs (X, z) . 

Theoreme 4.1.13. Pour £ > , la SGE E^(x) de la performance moyenne de I'algo- 
rithme glouton de couplage, quand I 'entree est choisie uniformement parmi les l-compo- 
santes, est telle que 



(4.66) 



y=T(x)b-i/(b-2)\ 

avec la fonction gi determinee, uniquement avec la condition ge(0) = , par une recurrence 

9ttv)Q + y) {l ^ mh - X) = UV). (4-67) 

Preuve. Nous obtenons, en differenciant (I4.65|) par rapport a z puis en y fixant alors 
z = 1 , une recurrence des E^( x) : 



K iE) E^(x) = ^-exp[b 



b\ 



1)1 



x,z, 



z=l 



soit 



T(x) 
bl 



■F®(x,l) 1 + 6 



^(-!)(X) 



T(x) 



T o T(x) + 



T(x) b d 
b\ ~d~z 



F®(x,z) 



z=l 



Soit 







az 



z=l 



Tlx) 



alors, 



j,»#(i)- (T ° t y ^(i) 



T{x) 
bl 



-F^(xA) 1 + 6 



■d x E(-V{x) 
T(x) 



T oT(x)) 



1 + 6 
1 + 6 



T(x) 
T(x) 



bl 

T{x) 1 



r o T(x) K W H t o T(x) + ±^_jr(<-i) ( X ) 



(4.68) 
(4.69) 

(4.70) 
(4.71) 

(4.72) 

(4.73) 
(4.74) 

(4.75) 



r o T(x) j VjQHt o T(x) + ^t^F^ o r o T(x)(4.76) 
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avec la fonction ^ telle que 

F^-V o T oT{x)=T{xf +1 F^- l \x). 

Nous obtenons alors 

^E^\x)- iTO b T _ { f\ x E^\x) = 

( x + b ^\x) T o \ kWhi o r(x) + 



+ 



T(x) 

roT(x) T p { t-i) OToT{x) 



b{b-l)T(x) 



Ainsi, la forme 



_T(x) a 



(4.77) 

(4.78) 
(4.79) 
(4.80) 

(4.81) 



(6-2)! 



permet le changement de variable 

r T( x ) b - 1 

' y 



(6-2)! 

^^(^^^(^(l + y)^ 6 - 1 ) 



(4.82) 



1 - 



sfc(y)(i + y 



-i) 



pour avoir 

v^vXi + v)** 6 "™ 6 - 1 ^ 



(l + by 



i-y 



g-i(y) \ t / \ , y s^- 1 ) 



avec 



T( X y# w {e) H e o T( 



(4.83) 



(y), (4.84) 



(4.85) 



r o T(x) 

Ce qui permet de conclure sur la recurrence des gg , done des E^ du theoreme en iden- 
tifiant 

<> 



^(y) = (i + ^ 



Proposition 4.1.14. La fonction go , telle que 

E^(x)=g (y) 

y=T(x) b - 1 /(b-2)\ 

admet un developpement en y = 1 gtii commence comme suit : 

. , (36-2)(l-l/2 b /( f> - 1 )) 
g0( ^ = 2% -y? ~ + - 



(4.87) 



(4.88) 
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Preuve. Pour I = , l'equation (|4.86p donne Jq telle que 



Nous avons 



w _ (1 + ^(l±^)_^M_ yWW + ^M (4 ,„) 

et a partir de (|4.85 j) , 

J ^ = ^W = W^y? + - (4 ' 91) 

En remarquant que 

5 _ 1 (y) = I((l + y) 1 /(^i)_l) , (4.92) 
nous avons aussi le developpement en y = 1 suivant : 

*(, l-» j (l + ^-l) = (1-,) +••• (4 - 93) 

Le premier terme du second membre de (I4.90P admet done un developpement en y = 1 
qui commence comme suit : 

(1 + (, l- y J (l +l/ )i/(»-i) )Jo(y) = 2(1 -y)3 + ' • • (494) 

Voyons maintenant comment commence le developpement en y = 1 du second terme, a 

savoir ^ (y) / (b 2 — b) , du second membre de (|4.90p . Si t = , dans les definitions (|4.77|) 
et (I4.72p nous obtenons : 

o r o r(a) = T(x)F ( -- 1 \x) (4.95) 

-bToT(x)$ x E(°\x). (4.96) 



2 = 1 



Dans T(x)-^F ( -°\x, z)\ , , les termes (des recombinaisons) qui contribuent au coefficient 
de (1/(1 — y)) de plus grande puissance dans le developpement de F^~ x \y) , en y = 1 , 
se deduisent des recombinaisons qui contribuent le plus avec les SGEs univariees. Le 
passage aux SGEs bivariees, puis 1' application de () | ont pour effet de changer 
exactement un terme $ x k Hj o T{x) en -d x k Ej{x) . Ici, pour I = , ces termes ou encore 
ces constructions proviennent de 

Tix) 1 1 (t ^-y ) 2 + ,) .or, 
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et sont 

1 „ , xx2 6(6-2) . A ^E^Xx) 
- (br o T(x)f + ^-^r o T(x)j T{x) 1 > . (4.98) 



Par le changement de variable 

T{x) b - 1 
(6-2)! 

nous avons 



V, (4.99) 



* m 2 Et-Q{x) (l-l/2^-D)(6-^ + ^ 



T(x) 6(1 -y)3 

Nous trouvons un developpement en y = 1 suivant : 

#(°>(y) _ (l-l/2 fc /( fe - 1 ))(6-l) 

6(6 -r ' - :t "' 

Ce qui, avec (|4.94|) , donnent 

My) = '"" — - + ... ( j..k>2> 

done 
et 

:/ni//i - — - + ... « 1 104) 





6(1 


-y) 3 


(36- 


-2)(1- 


!/ 2 6/(6-l)) 




26(1 - 




(36 


-2)(1- 


l/ 2 b/(6-l)) 




2 2 6(1 


-y? 


(36- 


-2)(1- 


1/2 6/(fe-l)) 




2 3 6((1 


-y) 2 



Proposition 4.1.15. L'equivalent asymptotique du coefficient [x n ] E^\x) de la SGE liee 
d la moyenne de la variable aleatoire Yn est : 

avec s , le nombre d'hyperaretes, et n = n{s) = s(b — 1) . 
Preuve. Notons E^ , la SGE lisse associee a la SGE #(°) . Mors, 

^ w ^o TW= ("-ff 1= y) + .... ^ 



et 



d m(f) _ (36-2)(l-l/2^))(6-l)r(Q 

dl E {t) ~ 226(1 -r(t)) 3 i + "" ( } 
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La formule d'inversion de Lagrange donne 



£ (0) (x) = (4.108) 
(b - 1)(» - 2)(1 - l/2V(f-D) J f r(t) + A e^df = 



2 3 6ivrn / V(l-T(i)) 3 "'/ *™ * 

(6 _ 1)(36 _ 2 )(l _ 1/2^-D) j f r(t) + A e"<«> dr(t) ^ 



2 3 6i7rn[(6-2)!] s (6- 1) / V(l-r(t)) 3 "7 r(t) s r(t) 
Nous en deduisons que 

[x n ]E^(x) = (4.111) 

(6-l)(36-2)(l-l/2 6 /(^D) 1 ( A exp^) df 2) 



dt. (4.113) 



2 3 6ivrn[(6-2)!] s / \(l-f) 3 

L'asymptotique de ce coefficient est done porte par 

(6-l)(36-2)(l-l/2 b /( fo - 1 )) f 1 exp(st) 



2 3 bi7rn[(b - 2)\} s J (1 - t) 3 t 



Cette expression a ete deja rencontree au facteur pres dans la preuve du lemme 13.3.31 . 
Par ce lemme, nous trouvons avec les valeurs £ = et m = 3 , l'equivalent asymptotique 
suivant : 

(6 _ i)(36 - 2)(1 - 1/2V(6-D) j 1 exp(*t)_ dt = 



2 3 6ivrn[(6-2)!] s J (1 - 1) 3 i 

(36-2)(l-l/2 f, /( f, - 1 ))(e/3) 3 / 2 e s / 1 



{l + 0(-^)}. (4.115) 



2%^[(b - 2)\} s 

Ainsi, nous concluons l'equivalent asymptotique du coefficient : 



Nous avons le coefficient [x n ] E<® (x) par cette proposition et l'asymptotique du coef- 
ficient [x n ] Hq o T(x) par la proposition 13 . 2 . ll . Nous trouvons alors par la definition de la 
SGE £(°)(s) : 

Theoreme 4.1.16. La moyenne En^ de la variable aleatoire Y^ correspondant a la 
taille de Vhypercouplage produit par I'algorithme glouton, quand V entree de Valgorithme 
est choisie uniformement parmi les hypercycles a n sommets, est telle que 

E f ~ , ( 4.w) 

avec s , tel que n = s(b — 1) . 
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Proposition 4.1.17. Pour £ > 1 , la fonction gi , telle que 

1 



T{x) 



,{y){l+yf^ 



y=T(x)»-l/(b-2)\ 



admette un developpement en y = 1 qui commence comme suit : 



9e ^ (3£ + 2)(1 - y) 3£ + 2 
ou bi admet la recurrence 

t (36-2)(l-l/2 fe /( fc - 1 )) 

bo ~ ¥b ' 

f (l-l/2 fe /( fc - 1 ))A, 
Cty = _»//,. — r 



+ ... , 



+ 



2^/C'- 1 )2 

(31 + 1)^2^)^ + 2 A^-x-^A^j 
2 £ /( fo - 1 )2 



(4.118) 

(4.119) 

(4.120) 
(4.121) 

(4.122) 



A,- etant la notation definie dans le theoreme \2.3.31 . 



Preuve. Dans la recurrence (|4.83p de la preuve du theoreme 14.1.131 . nous remarquons que 
pour determiner , il suffit d'integrer 



(l + by 



i-y 



(1+^)1/(6-1) )Ji(y) 



(l + y)(W-l)/(6-l) 



&(6-l)(l + y)(W-l)/(b-l) 



(4.123) 



Or, dans cette equation, le premier terme du second membre admet un developpement 
en y = 1 qui commence comme suit : 



{l + by 



l+(ft-l)fly 



{1 + y yi y (L) )My) (l-l/2 fe /( fe - 1 ))A^(6-l) 



2( 



car 



et 



(l + by 



(l+y)(^+6-l)/(6-l) 



3^+3 



1-2/ 



2^/(6-1)2(1 -y) 



5 _!(y) 1-1/2^-D 

_1)J . . r 



+ . 



(l+y)l/(6-l) 



(i-y) 



•My) = \~qtt~o + • • • 



(4.124) 
(4.125) 
(4.126) 



avec la notation A^ du theoreme 12.3.371 . Nous trouvons ensuite que le debut du develop- 
pement en y = 1 du second terme du second membre de (|4.123p est : 



6(6- 1)(1 +y)(*+6-l)/(&-l) " 

(6 _ i)M {(3^ + iJ&^-DAM) + 2 ^i-jbfP'W } 
2^/( fc - 1 )2(l -y) 3 ^+ 3 



(4.127) 
+ ... (4.128) 
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Comme dans le lemme 12.3.361 . si une composante £ subit la decomposition (suppression 
d' hyper aretes) suggeree par l'algorithme[3]produisant une famille, indexee par i , de com- 
posantes d'exces ji et ayant ki marquages (comptes avec l'ordre de multiplicity defini par 
le nombre d'hyperaretes supprimees auxquelles un sommet marque appartenait) alors, en 
supposant que (Ji, ki) ^ (— 1, 1) , le degre maximum de (1/(1 — y)) du developpement en 
y = 1 d'une SGE liee a la description de telles families, vaut au plus 

^(3j 4 + 2h) = (3£ + 3) - J2 k i ■ ( 4 - 129 ) 

i i 

Comme la SGE F , (|4.77|) , est liee a de telles families sauf a celle reduite a une composante 
d'exces £ et ayant un seul marquage, le degre est au plus, avec Yli k% = 2 , 

3^ + 3-2 + 2 = 3^ + 3, (4.130) 



3 ° 



ou le 2 rajoute est a cause d'un facteur en , d x k E^(x) remplacant un facteur en "d x k H, 
T(x) par la derivation M- dans (|4.72|) . Si b > 3 , nous distinguons deux types de families 
telles que Yli ki = 2 avec des ki d'ordre 1 (cette restriction est motivee parce que nous 
ne sommes interesses que par le coefficient du terme en (1/(1 — y)) avec la plus grande 
puissance) : 

- Si les deux marques appartiennent a une meme composante, alors les families sont 
liees a des constructions decrites dans 

C l { - t - 1) {x,z)= (4.131) 
T(x) + ^ [b j + - ^ j x (4.132) 

(x, z) . (4.133) 

- Si les deux marques appartiennent a deux composantes distinctes, alors les families 
sont liees a des constructions decrites dans 

C^- 1 \x,z) = 

T{x)£+1 . 1 f ^ x G^\x,z) b - 2 V | b^G^\x,zf- 




2 \ (fe-2)! I 2 (6-3)! 

e-i 

x ^ $xG {p) (x, z) x ■d x G^~ x ~ p>) (x, z) . (4.136) 

p=0 

D'un cote, l'application de () \ z= i a ()4.13ip nous permet de preciser les construc- 
tions contribuant au coefficient de (1/(1 — y)) portant la plus grande puissance. A partir 
de (|4.13ip , nous identifions alors les constructions decrites dans l'expression suivante 

~ (6r o T{x)) 2 + 6(6 ~ 2) r o T{x)\ T(x) e ~ 1 i!) x 2 E ( - e ~ 1 ' 1 (x) . (4.137) 
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Notons par bj(b — l) 2j /(3j + 2), le coefficient de (1/(1 — y)) portant la plus grande 
puissance de gj(y) dans le developpement de ce dernier, ce que nous avons ecrit 



9j(y) 



+ 



avcc 



(3j + 2)(l-y)3i+2 
(36-2)(l-l/2 b /( b - 1 )) 



(4.138) 
(4.139) 



Ces constructions, identifiers dans (|4.137p , contribuent done au coefficient par 

6(6 - 1) 2£+1 {U + l)k-i2 ie ' 1)/{b ' 1) • (4.140) 

De l'autre cote, l'application de () | 2= i a (j4.134H nous permet d'obtenir les autres 
constructions qui permettent de determiner ce coefficient de (1/(1 — y)) portant la plus 
grande puissance. Aussi, nous identifions les constructions decrites dans l'expression sui- 
vante : 



T(^(i(5roT(x)f + ^ 



toT(x) X 



-1 



x2^^EW(x)x^JJ w _ p or(x) 

p=0 



(4.141) 
(4.142) 



Ces constructions contribuent au coefficient par 

e-i 



26(6 -l^+^WA^ 1 ), 



(4.143) 



ou nous adoptons la notation A^ , du theoreme 12.3.371 . relative au coefficient du terme en 
(1 — t o T(x))~ 3i portant l'asymptotique du coefficient de Hi o T(x) : 



H e (t) 



i ( Xi{b-rr i 



if- \ U (1 - r(t)) 



3/' 



+ 



(4.144) 



avcc 





1 


A = 




2 


Ai = 


5 


8' 


\t = 


U 




~2 



(4.145) 



1 ^ 2 

+ - ^2 XpXi-i-p , pour 1 = 2, 3. 



P =i 



La contribution de F(y)/(b(b — 1)(1 + y)^ +h l ^/^ b ^) au coefficient est done 



2 £ /( 6 - 1 )2 



(4.146) 
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Nous obtenons alors par (|4.124|) et (|4.127j) la recurrence 
(1 _ 1/2 ^-i) )A 

b t ~ 2^=1)2 ( } 

(31 + l)j^2g-^ + 2 
+ • (4148) 



Ainsi, par exemple, nous trouvons : 
Exemple 4.1.18. Pour I = 1 , 



soit 



Exemple 4.1.21. Pour i = 4 , 



Proposition 4.1.22. La recurrence des be 



b\b - \) n 



s'ecrit encore 

2b t 2^) 
(1-1/2V(6-D) -^(6)3Af, 



<> 



( l- 1/2 ^)A 1 + 5So 



5 i2 i/(6-i) = ^ _ i/2&/(&-l))(7fc _ 4)^1 . (4.150) 
Exemple 4.1.19. Et pour £ = 2 , 

L 2 VM = (l^V^hh + 8£i2V(>-i)+2A 1 So 

^ 2 2' ' 

soit 

^ 2 2/(6-i) = ^ _ 1 / 2 V(t-i))(656 - 34)^| . (4.152) 
Exemple 4.1.20. Pour I = 3 , 

S 3 2 3 /( 6 " 1 ) = (1 - l/2 fc /(^ 1 ))(2936 - 144)^| . (4.153) 



g 42 4/(6-i) = (i _ i/2 6 /( & - 1 ))(70816 - 3314)^^ . (4.154) 



9t(v) = 7Wn I hr, > lW+ 2 + • • • , (4-155) 



(4.156) 
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avec \i , tel que 



An = - 



1 

2 ' 



X e = -A^_i(3^ - 1) + ^ ^2 X P X e-i- P , pour £ = 1,2. 



e-i 



(4.157) 



p=0 



et 



Mb) 



3b -2 
"36" ' 



(4.158) 



+^Z^ Xi - 1 -P X P 3^ ' pour £=1,2... 



p=0 



Proposition 4.1.23. Le coefficient de la SGE E^\x) qui code la moyenne de la variable 
aleatoire yJP admet V equivalent asymptotique suivant : 



x 



(b - l) 2i+1 (l - l/2 b /( b - 1 ))a e (b)3X e e s+M ^-^(s/(3£ + 3))( 3£ + 3 )/ 2 



4^i7r[(6-2)!] i 



x + 



(4.159) 
x (4.160) 

(4.161) 



Preuve. Comme la SGE lisse E^ correspondant a la SGE T(x) e E^ (x) est telle que 
F(t\ a oW/Vl +T ffflVH) {b-ir\l-l/2 b /^))a £ {b)3\, 



ct 



d_ . (6- !)*(! -l^- 1 ))^^, (ft-l)r(t) 
dt {> 2(1 - r(i)) 3 *+ 3 t 



(4.163) 



Par la formule d'inversion de Lagrange, l'asymptotique du coefficient [x n ] E^ (x) est 
porte par celui de 



i= (b- i)^+i(i -i^-D^^aA, 

2n 



h n+e 



r(t) 



(l-r(t))™ 



$(i) n ^ , (4.164) 



avec r(t) = r _ 7(6 - 2)! , $(t) = exp(r(t)/(6 - 1)) et m = 3£ + 3 . En notant n = n(a) = 
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s(b — 1) — t , nous obtenons le resultat par le lemme 13.3.31 car 

/ = (4.165) 

= 4nV55F[(6-2)!]- j 1 + °(^) j =( 4 - 166 ) 

( 6 - 1)^+1(1 _ l/2V(^i))g £ (6)3A,e- + ^-^i( S /(3£ + 3))(^+ 3 )/ 2 
= 4n^i5F[(6-2)!]- X (4 ' 167) 

x { 1 + °(^)} • ( 4 - 168 ) 

Par les asymptotiques des coefficients [x n ] E^\x) de la proposition precedente et [x n ] Hgo 
T(x) de la proposition 13.3.41 . nous trouvons : 

(£) , (£) 

Theoreme 4.1.24. La moyenne En de la variable aleatoire Yn correspondant a la 

taille de I'hypercouplage produit par I'algorithme glouton, quand V entree de I'algorithme 

est choisie uniformement parmi les hypercycles a n sommets est telle que 

avec s, le nombre d'hyperaretes quand n est le nombre de sommets, soit n = n(s) = 
s(b — 1) — I , et les definis dans I'enonce de la proposition \4-l-l 



4.2 Hypergraphes evoluants 

Dans cette section nous generalisons un resultat de [Hj sur les graphes evoluants. Plu- 
sieurs resultats sur les graphes aleatoires [5] seraient aussi accessibles pour hypergraphes 
par l'approche suggeree dans [14J . 



Algorithme 4 : Generation d'hypergraphe jusqu'a l'apparition du premier cycle 
Entrees : Entier n , le nombre des sommets etiquetes. 
Sorties : Foret d'hyperarbre non enracine et un cycle avec une hyperarete 
inter ieure marquee. 

debut 

Initialiser avec un hypergraphe a n sommets sans hyperarete. 
repeter 

I Ajouter une nouvelle hyperarete aleatoire 
jusqu'a une composante d'exces > a ete formee 

fin 

Marquer la derniere hyperarete ajoutee dans l'hypergraphe. 
retourner I 'hypergraphe 
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4.2.1 Definitions et notions 

Dans cette section, nous disposons de la SGE H^\(w,z) bivariee des hyperarbres : 

H i o T(w 1 ^ b ^ 1 h) 

n-iM= 1 ZZ-V ' (4170) 

avec 

- z : denote la variable liee au nombre de sommets de l'hypergraphe, 

- w : denote la variable liee au nombre de ses hyperaretes, 

- T : la SGE des hyperarbres enracines 

T(z) = zexp(^^) , (4.171) 



(6-1) 



H-i : la SGE "lisse" des hyperarbres 



H_ 1 (t) = t[l-^\ , (4.172) 



ou r(t) = . 
Definition 4.2.1. Denotons ^f n l'operateur 

*» F =E%f aveciV=(;), (4.173) 

m>l '"'Km) v 7 



et F(w, z) la SGE bivariee 



F(w,z)= V f m , n w mZ -. (4.174) 



m 

n,m>0 

Cet operateur \l/ n admet une formulation integrate suivante : 
Definition 4.2.2. l'operateur ty n est 

f°° 1 dt 

ou 

N = (f) (4.176) 



ei 

F n H = n\ [z n ] F(w, z) = ^i ^S^d * ■ (4-177) 

Inequivalence des deux definitions se justifie, comme souligne dans [Hj , par la j3 
fonction integrate : La substitution u = t/(l + t) transforme J °° t m ~ 1 dt/{\ + t) N+l en 



i i 
u m -\l- u) N - m d u = B(N + 1 - m, m) = . (4.178) 
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Definition 4.2.3. Un hypergraphe evoluant est un hypergraphe dont les hyperaretes 
apparaissent une par une. 

Ces hypergraphes evoluants sont a priori differents des hypergraphes ou chacune 
des hyperaretes existe avec une certaine probabilite p (modele binomial d'hypergraphe 
aleatoire etudier dans [10]) . 



4.2.2 L'attente moyenne de l'apparition du premier cycle 

Pour obtenir une expression de l'esperance du nombre d'hyperaretes ajoutees par 
l'algorithmeH], nous considerons la SGE bivariee liees aux structures rencontrees lors du 
deroulement de l'algorithme telles que le processus de creation d'une nouvelle hyperarete 
continu. Ces structures rencontrees ne sont autres que des forets d'hyperarbres, chacune 
de ces structures ayant m hyperaretes peut apparaitre en suivant ml historiques distinctes 
selon l'ordre de creation des m hyperaretes. Sachant qu'un hypergraphe ayant n sommets 
admet m hyperaretes ordonnees, la probabilite qu'il s'agit d'une foret d'hyperarbres est 

m lfm,n fm,n , . „„x 

avec N = m et avec f m>n le nombre de forets d'hyperarbres chacun ayant n sommets et 
m hyperaretes. Nous obtenons alors le nombre moyen d'hyperaretes ajoutees en sommant 
le terme (|4.179|) pour m > : 

E n = E ■ (4180) 

Notons F la SGE des forets d'hyperarbres 

F(w,z)=exp( ) {b _ 1} M= E f^w m ^, (4.181) 

\ ' / m>0,n>l 

alors l'application de l'operateur ^ n a F donne presque l'expression de la moyenne E n : 

*nF = £ %- , (4.182) 

m>l m U 

une expression plus proche de la moyenne E n est 

*»*« F = E ■ ( 4 - 183 ) 

m>l \m) 

Ainsi, nous voyons que la moyenne E n s'exprime avec l'operateur comme : 



E n = 1 + 



(4.184) 
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Notons alors que 

KF(w,z) = F(w,z) l wl/(6 _ 1)& / (4-185) 
T t w i/(b-i) y> {H- 1 oT(w 1 K b - 1 h)\ 

= U-HM eXP ( J<^> j ' (4186) 

Notons f n {w) , le coefficient [z n ] i9 w F(w, z) : 

r m 7V /(6_1) *) 6 /i?-ioT(^/(6-i)^)\ 
^ = ^ U^Dft! 6XP ( " J • ( 4 - 187 ^ 

Une application d'une version de la formule d'inversion de Lagrange 12.3.21 . readaptee a 
cette SGE bivariee avec une seconde variable relative au nombre d' hyper aretes, donne 

fnH = (4.188) 



l ^ {It (4 exp (*-^'))) ex p ( ntt ^)) = ( 4 - 189 ) 

j w ±. (^ e H-^,t)\ e n W ^-nHt) dti (4 . 190) 



H-i(w,t) =t I l-w^ ) . (4.191) 



2W7T 

avec r(t) = t^ 1 / (b - 2)\ , 



Nous utiliserons une autre formulation de la version de l'inversion de Lagrange qui ne 
comptabilise pas les hyperaretes en mettant la variable w du cote des forets, mais du 
cote de la structure centrale. Echanger la maniere de comptabiliser les hyperaretes pour 
les comptabiliser du cote de la structure centrale, se traduira par une multiplication par 
w n/(b-l) ^ n/(b — 1) etant le nombre d'hyperaretes dans une foret a n sommets pour 
garantir la presence d'au moins un cycle, et un changement de t en iu;" 1 /^ -1 ) . Nous 
obtenons alors une version de l'inversion de Lagrange suivant : 

fnH = (4.192) 

n Us \ \ > 

W»- 1 r.«_n/d / 1 t b ( =, , t A\ n m 



n 

n—1 



^ (h is>)) - (»^)) ■ <" M » 



n 

De maniere analytique, cela s'ecrit 

/.(«)- =£x (4.195) 

I (^exp (i?^, -!->)) exp („lfil - „M( t ))) dt . (4.196) 
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Nous avons 



t 



d ft 



dt \bl 



exp 



W b-1 



r(t) 



t 



6-1 



t 



(6-1)! 



+ (1 " r(t)) exp 



t 



r(t) 
b-1 



1 + — ^- (1 _ T (i)) ) (<XP 



iy 6-1 



r(t) 



r(t) 



(4.197) 
(4.198) 
(4.199) 
(4.200) 



Et f n (w) se met sous la forme suivante : 
/«(™) 



n-l 

10 



2ivrn(6 - 1) 



x ^r(t)(l + -4j-(l-r(t)))r' 



(4.201) 



^) e n ^i- nln{t) dt. (4.202) 



Pour evaluer ^ n , & w F qui s'ecrit sous la forme integrale suivante : 



1 , . dw 



o (l + w) N ny 'w(l + w) ' 



avec N 



et 



(4.203) 
(4.204) 



F n (w) = n\f n (w) , 

nous compensons ce facteur n! dans l'estimation de f n par le changement de variable 



W 



[wn) b 1 
dw^ dw 

w w , 

e b - 1 

n\f n (w) -» n\f n (- rj^j-) 

<wn)° 1 



(4.205) 



et nous obtenons 



(6- l)n! 



f f e& ~ 1 ' 



d u; . 



(4.206) 



Dans ce changement de variable (|4.205|) , nous avons 



J>-1 



,n-l 



At( / ) 
[wn)° 1 



2iirn n (b — l)w r 



r i ff (w,t)e n?t(t) dt. 



avec 



nwt 



g( w ,t)=r(t)(l + — (l-r(t)) 



(4.207) 
(4.208) 
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V«) = f^-K«) + ?(i-^). <«*») 

La derivee de cette fonction h est 

h'(t) = 1^-1 + ^(1- T (t)) = (1 - r(t))(^ - . (4.210) 
La derivee seconde de la fonction /i est 

h"{t) = -{b - - ~) + (1 - r(t))(i) . (4.211) 

Nous distinguons alors deux cas : 

1. Si w < [(b_ 2 )!p/(b— i) ' a l° rs la derivee seconde au point col to = [Q> — 2)!] 1 /( b ~ 1 ) est 
positive 

h "^ = -'"-') ,(„ - 2) !]^-i) ( 7 - j(i - 2)V"-'» <4 ' 212) 

Alors nous faisons passer le contour d'integration verticalement par le point col to ■ 

2. Si u> > Tn^JMT/{b=Tj ; alors la derivee seconde au point col ti = e/w est positive 



e^ 1 s w 2 



Alors nous faisons passer le contour d'integration verticalement par le point col t\ . 
Dans chaque cas, nous utilisons le contour d'integration suivant selon tj : 

f 7l : t = tj + iv/n , v |G [Sntj, 3ni,] . . 

' ' \ 70 : t = tj + 3tj exp(ia) , a |e [vr/2, 3tt/2] . \ ■ J 

Notons 7i , une portion du chemin 71 : 

71 :t = tj + iv/n, v |e [— « s , k s ] , (4.215) 

avec k s , un nombre positif qui sera precise plus tard. Soit alors la valeur Ix de Pintegrale 
()4.207p restreinte a cette portion 71 : 

-j- — r / g(w,t)e nh ®dt. (4.216) 



Ce qui, si Pintegrande est approchee au point col, donne 

7i = 2^(b-i)Vi A exp - ^ )2/2 + °(^ w c*i)c - ^) 4 ) dt ' ( 421? ) 
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soit, avec la variable d'integration v 



Soit, encore par le changement de variable 



r 1 / /i"fe) r 



(4.219) 



et en notant 

C w , n = y/h»(tj)/(2n) , (4.220) 

h = (4.221) 

1 fl(ll! t ■\p nh ( t i) rK s C w , n r 2 | Q(- ft(4)(t J ) r 4 \ 

9{w,t,)e I ~ 1+U[ ^ T - W ™ ) dr= (4.222) 



vm™(6 - l)#-V2n/i"(^) J- Ks C w , n 

e n - l g{w,t j )e nh ^ r**C„, n / frC 4 )^ ■■' 



r " % " r i i+o( ^^ dr=(4223) 

(4.224) 



vrn«(6- l)#-V2n/i"(i j ) i- Ks c m ,„ V K h"(tj) 2 n 



irn n (b - l)w n - 1 y / 2nh"(t j ) 
Jf°"^r + 0(^M(^,). ,4,25, 



Alors l'integrale Ji devient 



h = yy , y (4.226 

7rn n (6 - \)w n ~ x ^j2nh!'{tj) 



x 



e~ r dr + Oih^yJhf'itj)^) . (4.227) 



En prenant k s = n 2//3 , nous obtenons 



(4.228) 



Sur le contour restant, la contribution etant exponentiellement petit, nous obtenons que 
I\ est l'ordre asymptotiquc dc (|4,207p : si w < e/[(b — 2)!] 1 /( b ~ 1 ) , alors on prendra j = 
soit le point col to = [(b — 2)!] 1 /( b_1 ) , sinon on prendra j = 1 soit le point col t\ = e/w . 
Et l'ordre asymptotique de (|4.207p est 
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Ann d'estimer asymptotiquement l'esperance, nous portons le terme asymptotique prin- 
cipal de l'equation precedente dans (|4.206p et nous avons alors a considerer l'integrale 
ainsi obtenue suivante : 

n\e n ~ l f°° g(w,tj) e nh{t 3 )-n\n{w) 



f 00 a(w,ti) e ntHtj)-nta{v>) 

/ ^ V h \ 3) dw, (4.230 



avcc 



n 



no-m^ ^<-^w^ (4 . 231) 

[ e/w siw> [(6 _ 2)l] i /(6 -i) , 

g(w, t) = r(t) (l + ^ (1 - r{t)) \ , (4.232) 

h{t) = ^ - ln(t) + ^(1 - , (4.233) 
la derivee seconde de cette dernier e etant 

h"(t) = -{b - - \) + (1 - r(i))(^) . (4.234) 

Pour evaluer cette integrate (|4.23U|) . nous distinguerons les deux intervalles ]0, Ao[ et 
]Ao, oof avec Ao = e/((6 — 2)!) 1 /( b_1 ) . Soit alors I a la valeur de l'integrale sur ]0, Ao[ : 

rl \ f ,n—l f\o i p nh(tj)—nln(w) 
I a = , / ^^ dw, (4.235) 



ou 



et 



w x 1 i , n (b-l)tjiv 6-2 , , (b-l)w 

m = 5^ - infe) + = + ln(Ao) + L_J_ (4.236) 



2 

1 wt,- v „ ^ / Aq \ ^ w 



= - l)ra (1 - — ) = ( b - 1) - (1 - 7") ■ ( 423T ) 
tj e \ e / Ao 

Cette derniere admet une racine en w = Ao ■ Par les equivalents suivants : 

iV + l-ra 6 - 1 ^ (4.238) 

et 

« 1 + (^""'>" /H ~ eS( " r ' < 4239 » 



4.2. HYPERGRAPHES EVOLUANTS 



93 



nous obtenons 



1 f X ° e'^iwt 1 +^(ij)-n\n(w) 



n\e n " / " u e 



h ~ ' ; / ; dw (4.240) 

n!e«-M(6-2)!]^TAo x 



n n (6- l)5Vvr2n 

"Ao e -^(^) i_1 +n¥^-»fe!-nln(^) d ^ 



n!e™- 1 [(6- 2) !]("-!) A f 1 e ^-D 



(4.242) 



n n (6- l)5Vvr2n Jo y/l-w 

"1 R nh(w) 



-dw (4.243) 



avec /i , definie par 



du), (4.244) 



! /1\ 6_1 6-1 6-2 



" W = -Ki^T)U; + — •" — -'■><-), (4.245) 
qui s'annule pour w = 1 . La derivee de cette fonction est 

6 \u; J b w 
qui s'annule pour la valeur de u> = 1 . La derivee seconde de h est 

/ 1 \ 6+1 1 

h"(w) = -(-) + — , (4.247) 

qui s'annule pour la valeur w = 1 done le point col est d'ordre 2 . La derivee troisieme de 
h est 

n &+2 2 



fe'» = (6 + l) - - — , (4.248) 

done /?/"(l) > (1 est un point d'inflexion de h qui est croissante) . Notons alors l'integrale 
I' a , la restriction de de l'integrale I a dans un voisinage de 1 , suivante : 



= !LL / f ^= dw (4.249) 

n n (b- l)2Vvf2n Ji-e n \/l - w 

r) \ p n re n -n((b-l)w 3 /6+0(w 4 )) 

= — x / = dw (4.250) 

n ra (6- l)2 V / vr2n V w 

= — x / = (l + 0(nw 4 )))dw. (4.251) 

n n (b- 1)5^2^70 Vw ' 
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Ceci, par le changement de variable 



w — 
dw 



(n(6-l)^V3 (4 252) 



donne 

„| p n fi i / -e„(n(b-l)/6)V3 p _ r 3 / ^ 4 N 

.253) 



n n (b- l)an6^2^ Jo V r V n 1 ^ y 



En prenant e n = n 1//4 nous obtenons 



nV6 

Nous en deduisons que la valeur asymptotique de I a est : 



h = 0{-^). (4.255) 



n 



La moyenne recherchee est alors determinee par revaluation de la valeur de l'integrale 

(|4.23Up restreinte a l'intervalle ]Ao, oo[ . Notons I a , cette integrale sur cet intervalle : 

n\e n ~ l f°° g(w,e/w) e nfofe)-nl n (w) 

I a = ; / r— i , — aw, (4.256) 

ou 

= ^4^2Y. (* + T - ^W^y)) ■ (4257) 

M^) = ^y+lnH, (4.258) 

et 

fe// fe) = +( 1 - ^-i ( V-2)! )( ^ ) - (4259) 

Cette derniere admet une racine en u; = A = e /((6 - 2)!) 1 /^ 1 ) . Par les equivalents 
suivants 

iV + l^n^ 1 ^ (4.260) 

et 

((l + T^rr^-efi^r, (4.261) 
(wn)° 1 



nous obtenons 



n , e n-l poo 



I a ~ / !ha e -g(f) M ^fe)-^W d , (4.262) 

n!e n Z" 00 g(w,e/w) n ( e ) b ~ l \ n ^ ? " 

~ == / 4= p +n ^ EirT ^ . 4.263 

ra n Vvr2n An /i e^ 1 w 
1 ' 1 w&-l(6-2)! 
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Ceci, par le changement de variable 



w 



{ W {b-2)\y/^- 1 ) 

dw dw 



w (6 — l)w 

g(w, e/w) -> w (l + ^ (1 - wfj , 



donne 

L 

Soit, 



me 



n" 



1 i i n_n 
b b 



W 



— w n n (b — l)Vvr2n \ 36 



n!e n / 66 + 2n 



(4.264) 



(4.265) 



Theoreme 4.2.4. Le nombre moyen d'hyperaretes au moment de I 'apparition du premier 
cycle dans un hypergraphe evoluant a n sommet est asymptotiquement 



2n 



a 36(6 - 1) 



(4.266) 
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Chapitre 5 

Annexe 



5.1 Preuves des deux identites combinatoires 

Lemme 5.1.1. Pour j, a G IN* (done j + a > 0) , 

i=0 v 7 i=0 v J 7 



Preuve. 



= a-»r(vy(.^r-|(-1^ + ge + ;:r , )«- 



t=J X 7 !=() X J 



II suffit de montrer 



htt}^ - £(' + ;_T>-^ 
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Comme 



j + a — i — 1\ x - 1 1 



ct 



j + A' + a-fc-lA fi + a-k-2\ ( 'j + a - k - 2\ f j + a-k-2\ 

(5.13) 

k 

j + a\ fj + a-k-l\ ^fj + a-k + r-l\ 

\ i-l i + r ) <514) 

nous obtenons l'identite. 
Lemme 5.1.2. Si a — j > alors 

e^ = jz( a ~^~ l )(i-ef-\ (5.15) 

i=o ^ J ' 

oil si k £ ~\N et t £ 2L alors 

't\ - { 1 \ _*(*-!)••• (*- k + 1) 



fey u-fey fe! 



Preuve. 



r=0 



(5.16) 



i=0 ^ ' 
f r -j - 1 
r 



= £ r ) o- - °y (5.i8) 

= Efi)(-l) r (l-e) r (5-19) 



= £ u-ira-^r (5.20) 

r=0 

= (1-(1 -*))''• (5.21) 

<> 
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LALGORITHME GLOUTON D'HYPERCOUPLAGE 
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5.2 L'hypergraphe utilise pour le deroulement de l'algo- 
rithme glouton d'hypercouplage 

#les 55 hyperaretes de l'hypergraphe 4-uniforme 

27 29 34 10 
22 12 32 30 
29 18 4 27 

22 15 3 18 
18 3 13 2 
4 23 8 33 
26 13 28 3 
32 17 31 1 
11 7 32 21 
6 26 4 22 

18 26 3 33 
21 2 12 23 

6 12 8 4 

4 22 10 19 

19 27 15 20 
15 9 26 
32 30 34 

20 19 2 
17 24 33 19 

13 33 26 25 

7 29 6 17 

14 2 9 6 
34 25 18 8 

28 27 3 26 
28 14 23 12 
17 23 6 21 
26 31 19 29 

21 3 13 16 

9 31 7 11 
17 13 11 19 
32 31 13 24 
4 17 14 

4 6 31 27 

13 14 25 16 

10 9 22 13 
2 34 31 
9 23 14 22 

11 33 31 7 

23 21 7 4 
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15 34 8 22 

30 12 22 
25 28 9 12 

15 30 2 27 
32 1 24 6 

24 15 6 

25 2 8 14 
27 18 28 

26 4 15 34 
4 33 7 32 
19 24 15 33 

16 12 34 17 

31 1 9 2 
13 16 33 
26 17 21 10 
29 25 2 5 



Bibliographie 



[1] Cayley A. A theorem on trees. Quart. J. Math. Oxford Ser., 23 :376-378, 1889. 

[2] Joyan A. Une theorie combinatoire des series formelles. Adv. in Math, 42 :l-82, 
1981. 

[3] Renyi A. On connected graphs i. Publ. Math. Inst. Hungarian Acad. Sci., 4 :385-388, 
1959. 

[4] Voblyi V. A. Wright and stepanov-wright coefficients. Mathematical Notes, 42 :969- 
974, 1987. 

[5] Bollobas B. Random Graphs. Academic Press, London, 1985. 

[6] Canfield E. R. Bender E. A. and McKay B. D. The asymptotic number of labeled 
connected graphs with a give number of vertices and edges. Random Structures and 
Algorithms, 1 :127-169, 1990. 

[7] Canfield E. R. Bender E. A. and McKay B. D. Asymptotic properties of labeled 
connected graphs. Random Structures and Algorithms, 3 :183-202, 1992. 

[8] Norman Bleistein and Richard A.Handelsman. Asymptotic Expansions of Integrals. 
Dover, 1986. 

[9] Berge C. Graphs and hypergraphs. North-Holland Mathematical Library, 1976. 

[10] Amin Coja-Oghlan, Cristopher Moore, and Vishal Sanwalani. Counting connected 
graphs and hypergraphs via the probabilistic method. 

[11] Martin Dyer, Alan Frieze, and Boris Pittel. The average performance of the greedy 
matching algorithm. The Annals of Applied Probability, 3 :526-552, 1993. 

[12] Knuth D. E. The Art Of Computer Programming. 2nd Edition, Addition- Wesley, 
1973. 

[13] Philippe Flajolet and Andrew M. Odlyzko. Random mapping statistics. 

[14] Knuth D. E. Flajolet P. and Pittel B. The first cycles in an evolving graph. Discrete 
Math, 75 :167-215, 1989. 

[15] De Bruijn N. G. Asymptotic Methods in Analysis. Dover, New- York, 1982. 

[16] Labelle G. Une nouvelle demonstration combinatoire des formules d'inversion de 
lagrange. Adv. in Math, 42 :217-247, 1981. 

[17] Selivanov B. I. Perechislenie odnorodnykh hipergrafov c protstoi ciklovoi strukturoi. 
Kombinatornyi Analiz 2, pages 60-67, 1972. 



101 



102 



BIBLIO GRAPHIE 



[18] Luczak T. Janson S., Knuth D. E. and Pittel B. The birth of the giant component. 
Random Structures and Algorithms, 4 :233-358, 1993. 

[19] Karohski M. and Luczak T. The number of sparsely edges connected uniform hy- 
pergraphs. Discrete Math, 171 :153-168, 1997. 

[20] Karohski M. and Luczak T. Random hypergraphs, 2000. 

[21] Karohski M. and Luczak T. The phase transition in a random hypergraph. J. 
Comput. Appl. Math., 142 : 125-135, 2002. 

[22] Wright E. M. The number of connected sparsely edged graphs. Journal of Graph 
Theory, 1 :317-330, 1977. 

[23] Wright E. M. The number of connected sparsely edged graphs, ii. smooth graphs 
and blocks. Journal of Graph Theory, 2 :299-305, 1978. 

[24] Wright E. M. The number of connected sparsely edged graphs, hi. asymptotic results. 
Journal of Graph Theory, 4 :393-407, 1980. 

[25] Bagaev G. N. Random graphs with degree of connectedness equal 2. Discrete Ana- 
lysis, 22 :3-14, 1973. 

[26] Flajolet P. and Sedgewick R. Analytic combinatorics. To 

appear chapters are available as Inria research reports. See 
http :/ /algo. inria.fr/flajolet/Publications/books. html. 

[27] Andriamampianina T. and Ravelomanana V. Enumeration of connected uniform 
hypergraphs, 2005. 

[28] Ravelomanana Vlady. Graphes multicycliques etiquetes : Aspects combinatoires et 
probabilistes, 2001. 

[29] Moon J. W. Various proofs of cayley's formula for counting trees. A seminar on 
graph theory, pages 70-78, 1967. 

[30] Herbert S. Wilf. generatingfunctionology. Academic Press, 1990. 



Hypergraphes aleatoires et algorithmiques 



Resume : Les hypergraphes sont des structures decomposables ou descriptibles done 
peuvent etre enumeres recur sivement. Ici, avec les fonctions generatrices exponentielles, 
nous obtenons des resultats d'enumerations exactes et asymptotiques des hypergraphes 
connexes a nombre de sommets et a nombre d'hyperaretes donnes. Dans un cadre combi- 
natoire, par un raisonnement d'inclusion exclusion, nous aboutissons a un encadrement 
des nombres des composantes d'hypergraphes : e'est une generalisation de l'encadrement 
de Wright pour les graphes. Pour obtenir les resultats asymptotiques, la methode du 
point col permet, en passant par l'analyse complexe, d'obtenir des demonstrations qui 
sont au final tres lisibles grace a Putilisation des fonctions generatrices. Soulignons que 
nous avons ainsi caracterise : 

- les composantes a nombre de sommets et a nombre d'hyperaretes donnes par rap- 
port a la taille moyenne d'un hypercouplage aleatoire de ces structures, 

- les hypergraphes aleatoires (evoluant hyperarete par hyperarete) par rapport au 
nombre moyen d'hyperaretes pour l'apparition du premier cycle. 

Cette these laisse envisager la possibilite de mieux connaitre les phenomenes de seuil avec 
des hypergraphes, ceci en s'inspirant des lignes de preuves qui s'y trouvent. 

Mots-cles : Hypergraphes uniformes, enumeration exacte, enumeration asymptotique, 
hypcrcouplages aleatoires, hypergraphes aleatoires, hypergraphes evoluants, analyse com- 
binatoire, fonctions generatrices, formule d'inversion de Lagrange, encadrement par in- 
clusion exclusion, methode du point col. 



Random hypergraphs and algorithmics 

Abstract : Hypergraphs are structures that can be decomposed or described ; in other 
words they are recursively countable. Here, we get exact and asymptotic enumeration 
results on hypergraphs by means of exponential generating functions. The number of hy- 
pergraph components is bounded, as a generalisation of Wright inequalities for graphs : 
the proof is a combinatorial understanding of the structure by inclusion exclusion. Asymp- 
totic results are obtained, proofs are at the end very easy to read thanks to generating 
functions, through complex analysis by saddle point method. We characterized : 

- the components with a given number of vertices and of hyperedges by the expected 
size of a random hypermatching in these structures. 

- the random hypergraphs (evolving hyperedge by hyperedge) according to the ex- 
pected number of hyperedges when the first cycle appears in the evolving structure. 

This work is an open road to further works on random hypergraphs such as threshold 
phenomenon, for which tools used here seem to be sufficient at first sight. 
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